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This review describes the properties of hadronic phases of dense matter in compact stars. The 
theory is developed within the method of real-time Green's functions and is applied to study of 
baryonic matter at and above the saturation density. The non-relativistic and covariant theories 
based on continuum Green's functions and the T-matrix and related approximations to the self- 
energies are reviewed. The effects of symmetry energy, onset of hyperons and meson condensation 
on the properties of stellar configurations are demonstrated on specific examples. Neutrino in- 
teractions with baryonic matter are introduced within a kinetic theory. We concentrate on the 
classification, analysis and first principle derivation of neutrino radiation processes from unpaired 
and superfluid hadronic phases. We then demonstrate how neutrino radiation rates from various 
microscopic processes affect the macroscopic cooling of neutron stars and how the observed X-ray 
fluxes from pulsars constrain the properties of dense hadronic matter. 
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1 Introduction 

Neutron stars are born in a gravitational collapse of luminous stars whose core mass exceeds the Chan- 
drasekhar limit for a self-gravitating body supported by degeneracy pressure of electron gas |T]. Being 
the densest observable bodies in our universe they open a window on the physics of matter under ex- 
treme conditions of high densities, pressures and strong electromagnetic and gravitational fields. Most 
of the known pulsars are isolated objects which emit radio-waves at frequencies 10*^ — 10^° Hz, which 
are pulsed at the rotation frequency of the star. Young objects, like the pulsar in the Crab nebula, are 
also observed through X-rays that are emitted from their surface as the star radiates away its thermal 
energy. Relativistic magnetospheres of young pulsars emit detectable non-thermal optical and gamma 
radiation; they could be sources of high-energy elementary particles. Neutron stars in the binaries are 
powered by the energy of matter accreted from companion star. 

The radio observations of pulsars stretch back to 1967 when the first pulsar was discovered. Since 
then, the observational pulsar astronomy has been extremely important to the fundamental physics and 
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astrophysics, which is evidenced by two Nobel awards, one for the discovery of pulsars (Hewish 1974) [2], 
the other for the discovery of the first neutron star - neutron star binary pulsar (Hulse and Taylor, 
1993) jS] whose orbital decay confirmed the gravitational radiation in full agreement with Einstein's 
General Theory of Relativity. The measurements of neutron star (NS) masses in binaries provide one 
of the most stringent constraint on properties of the superdense matter. The timing observations of 
the millisecond pulsars set an upper limit on the angular momentum that can be accommodated by a 
stable NS - a limit that can potentially constrains the properties of dense matter. Noise and rotational 
anomalies that are superimposed on the otherwise highly stable rotation of these objects provide a clue 
to the superfluid interiors of these stars. 

The orbiting X-ray satellites allow astronomers to explore the properties of the NS's surface, in par- 
ticular its composition through the spectrum of thermal radiation and the transients like thermonuclear 
burning of accreted material on the surface of an accreting NS. The currently operating Chandra X-ray 
satellite, Newton X-ray Multi Mirror Mission (XMM), Rossi X-ray timing explorer (RXTE) continue 
to provide new insight into compact objects, their environment and thermal histories. Confronting 
the theoretical models of NS cooling with the X-ray observations constrains the properties of dense 
hadronic matter, its elementary particle content and its condensed matter aspects such as superfluidity 
and superconductivity. 

The currently operating gravitational wave observatories VIRGO, LIGO, GEO and TAMA are 
expected to detect gravity waves from various compact objects, the NS-NS binaries being one of the 
most important targets of their search. Higher sensitivity will be achieved by the future space-based 
observatory LISA, currently under construction. The global oscillations of isolated NS and NS in binaries 
with compact objects (in particular with NSs and black holes) are believed to be important sources of 
gravity waves, which will have the potential to shed light on the internal structure and composition of 
a NS. 

The theory of neutron stars has its roots in the 1930s when it was realized that self-gravitating matter 
can support itself against gravitational collapse by the degeneracy pressure of fermions (electrons in the 
case of white dwarfs and neutrons and heavier baryons in the case of neutron stars). The underlying 
mechanism is the Pauli exclusion principle. Thus, unlike the ordinary stars, which are stabilized by 
their thermal pressure, neutron stars owe their very existence to the quantum nature of matter. The 
theory of neutron stars has been rapidly developing during the past four decades since the discovery of 
pulsars. The progress in this field was driven by different factors: the studies of elementary particles 
and their strong and weak interactions at terrestrial accelerators and the parallel developments in 
the fundamental theory of matter deeply affected our current understanding of neutron stars. The 
concepts of condensed matter, such as superfluidity and superconductivity, play a fundamental role in 
the dynamical manifestation of pulsars, their cooling and transport properties. Another factor is the 
steady increase of computational capabilities. 

This review gives a survey of the many-body theory of dense matter in NS. It concentrates on the 
uniform phases and hadronic degrees of freedom which cover the density range Po < p < 3po, where po 
is the nuclear saturation density. The theory is developed within the framework of continuum Green's 
function technique at finite temperature and density jH IHl Ej • Such an approach allows for a certain 
degree of coherence of the presentation. Subsection 12 . 1 Ul at the end of Chapter 2 gives a brief summary 
of methods and results omitted in the discussion. It is virtually impossible to cover all the aspects 
of the NS theory even when restricting to a certain domain of densities and degrees of freedom. The 
topics included in this review are not surprisingly aligned with the research interests of the author. It 
should be noted that there are very good textbooks |H1 E], monographs [TTl [T21 CHI and recent lecture 
notes J3l that cover the basics of compact star theory much more completely than it is done in this 
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review. Furthermore there are several recent comprehensive reviews that cover more speciahzed topics 
of NS theory such as the role of strangeness in compact stars [14^, cooling of neutron stars f^l ITBj . 
neutrino propagation jTT] , phases of quark matter at high densities [IHl CSl 120] and the equations of 
state of hadronic matter j5Tl 1^ 1^. 

The remainder of this introduction gives a brief overview of phases of dense matter in neutron 
stars. Chapter 2 starts with an introduction to the finite temperature non-equilibrium Green's functions 
theory. Further we give a detailed account of the T-matrix theory in the background medium and discuss 
the precursor effects in the vicinity of critical temperature of superfluid phase transition. Extensions of 
the T-matrix theory to describe the superfluid phases and three-body correlations are presented. A virial 
equation of state is derived, which includes the second and the third virial coefficients for degenerate 
Fermi-systems. The formal part of this Chapter includes further a discussion of the Brueckner-Bethe- 
Goldstone theory of nuclear matter, the covariant version of the T-matrix theory (known as the Dirac- 
Brueckner theory). We further discuss isospin asymmetric nuclear matter, onset of hyperons and meson 
condensates, and the conditions of charge neutrality and /^-equilibrium that are maintained in compact 
stars. We close the chapter by a discussion of stellar models constructed from representative equations 
of state. 

Chapter 3 discusses the weak interactions in dense matter within the covariant extension of the 
real-time Green's functions formalism described in Chapter 2. We discuss the derivation of neutrino 
emissivities from quantum transport equations for neutrinos and the dominant processes that contribute 
to the neutrino cooling rates of hadronic matter. Special attention is paid to processes that occur in 
the superfluid phases of neutron stars. 

Chapter 4 is a short overview of cooling simulations of compact stars. It contains several example 
simulations with an emphasis on the processes that were discussed in Chapter 3. 

1.1 A brief overview of neutron star structure 

We now adopt a top to bottom approach and review briefly the sequence (s) of the phases of matter 
in neutron stars as the density is increased. A schematic picture of the interior of a M = I.4M0 
mass neutron star is shown in Fig. ^ The low-density region of the star is a highly compressed fully 
ionized matter at densities about p = 10^ g cm^ composed of electrons and ions of ^^Fe. This phase 
could be covered by several cm thick 'blanket' material composed of H, He, and other light elements, 
their ions and/or molecules. The composition of the surface material is an important ingredient of the 
photon spectrum of the radiation, which is used to infer the surface photon luminosities of NS. Charge 
neutrality and equilibrium with respect to weak processes imply that the matter becomes neutron rich 
as the density is increased. In the density range 10^ ^ P ^ 10^^ g cm~^ a typical sequence of nuclei that 
are stable in the ground state is ^^Ni, ^^Kr, ^^Se, ^^Ge, ^°Zn, ^^^Mo, ^^^Zr, ^^°Sr and their neutron rich 
isotopes. The matter is solid below the melting temperature T^ ~ 10^ K, the electron wave functions 
are periodic Bloch states, and the elementary excitations are the electron quasiparticles and phonons. 
The lattice may also contain impurities, i. e. nuclei with mass numbers different from the predicted 
stable nucleus, as the time-scales for relaxation to the absolute ground state via weak interactions could 
be very large. The transport properties of the highly compressed solid (CS) are fundamental to the 
understanding of the way the thermal energy is transported from the core to the surface and the way 
the magnetic fields evolve in time. 

Above the density p ~ 4 x 10^^ g cm~^ not all the neutrons can be bound into clusters, and those 
which are free to form a continuum of states fill a Fermi-surface characterized by a positive chemical 
potential. Thus, the "inner crust" is a CS featuring a neutron fiuid. The sequence of the nuclei in the 
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R = 10 [km] 




M=1.4Mq logp = l [g/cm^] 

Figure 1: Schematic cross section of a M = IAMq mass neutron star. The different stable 
phases and their constituents (using the standard notations) are shown as a function of the 
radius in the upper half of the diagram and as a function of the logarithm of the density p 
in the lower half. Note that the both scales are strongly expanded in the low-density/large- 
radius domain. For details see the text. 



inner crust are neutron rich isotopes of Zr and Sn, with the number of protons Z = AO and 50 respectively 
and mass number in the range 100 < A < 1500. Below the critical temperature ~ 10^° K, which 
corresponds to 1 MeV (1 MeV = 1.602 x 10^° K) neutrons in the continuum undergo a phase transition 
to the superfluid state. At low temperatures, the electron quasiparticles and the lattice phonons are the 
relevant degrees of freedom which control the thermal and magnetic properties of the matter. At non- 
zero temperatures neutron excitations out of condensate can play an important role in mass transport 
and weak neutral current processes. 

At about half of the nuclear saturation density, po = 2.8 x 10^^ g cm~^ the clusters merge into 
continuum leaving behind a uniform fluid of neutrons, protons and electrons. The uniform neutron (n), 
proton (p) and electron (e) and possibly muon (/i) phase extends up to densities of a few po! the p, e 
and n abundances are in the range 5-10%. The many-body theory which determines the energy density 
of matter in this density range (as well as at higher densities) is crucial for the structure of the neutron 
stars, since most of the mass of the star resides above the nuclear saturation density. The neutrons 
and protons condense in superfluid and superconducting states below critical temperatures of the order 
10^ K. Because of their low density the protons pair in the relative ^5*0 state; neutron Fermi-energies lie 
in the energy range where attractive interaction between neutrons is in the F2 tensor spin-triplet 

channel. The relevant quasiparticle excitations of the npep-phase are the electrons and muons at low 
temperatures; at moderate temperatures the neutron and proton excitations out of the condensate can 
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Figure 2: The Schwinger (left) and Keldysh (right) contours; the axis corresponds to the real 
time, to is the observation time. The time flows from left to the right on the upper branch 
and from right to the left on the lower branch of the contour. 



be important. 

The actual state of matter above the densities 2 — 3 x po is unknown. The various possible phases 
are shown in Fig^ At a given density the largest energy scale for charge neutral and charged particles 
are the Fermi-energies of neutrons and electrons, respectively. Once these scales become of the order 
of the rest mass of strangeness carrying heavy baryons, the S^'°, A, hyperons nucleate in matter. 
Their abundances are again controlled by the equilibrium with respect to weak interactions and charge 
neutrality. Since the densities reached in the center of a massive NS are about lOxpo, it is likely that 
the critical deconfinment density at which the baryons lose their identity and disintegrate into up (m), 
down {d) , and possibly strange (s) quarks is reached inside massive compact stars. The critical density 
for the deconfinement transition cannot be calculated reliably since it lies in a range where quantum 
chromodynamics (QCD) is non-perturbative. NS phenomenology is potentially useful for testing the 
conjecture of high-density cold quark matter in compact stars. 

At densities p > 2po Bose-Einstein condensates (BEG) of pions (vr), kaons (K), and heavier mesons 
can arise under favorable assumptions about the effective meson-nucleon and nucleon-nucleon interac- 
tions in matter. For example, the pion condensation arises because of the instability of the particle-hole 
nucleonic excitations in the medium with quantum numbers of pions. This instability depends on the 
details of the nuclear interaction in the particle-hole channel and is uncertain. There are very distinct 
signatures associated with the pion and kaon BEG in the physics of NS featuring such a condensate, 
which includes softening of the equation of state and fast neutrino cooling. 

2 The nuclear many-body problem 

2.1 Real-time Green's functions 

Consider a non-relativistic Fermi system interacting via two-body forces. The Hamiltonian in the second 
quantized form is 



where are the Heisenberg field operators, x = {r,t) is the space-time four vector, a stands 

for the internal degrees of freedom (spin, isospin, etc.), m is the fermion mass and V{x,x') is the 
interaction potential, which we assume to be local in time V{x,x') = V{r,r')6{t — t'). The creation 
and annihilation operators obey the fermionic anti-commutation rules, {ilJa'{x'),4'l{x)} = S{x — x') and 
{ipai{x'),ilJa{x)} = 0, and their equation of motion is given by idtip{x) = [iIj{x),7{\ , where [, ] and {, } 
stand for a commutator and anti-commutator (here and below h = c = 1). The fundamental object of 
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the theory is the path-ordered correlation function 

G(l,l') = -^W(l)^ni')), (2) 

where the path-ordering operator arranges the fields along the contour such that that time arguments 
of the operators increase from left to right as one moves along the contour shown in Fig. |21 (here and 
below the boldface characters stand for functions that are ordered on the contour). The unitary time 
evolution operator propagates the fermionic wave function according to ^'(^2) = S{ti,t2)ipiti) along the 
Schwinger contour from —00 — >■ to ~^ +00, where to is the observation time (Fig. |2)). The Keldysh 
contour is obtained from the above one by inserting a piece that propagates from to —>■ +00 and 
back. The Keldysh formalism is based on a minimal extension of the usual diagrammatic rules, where 
depending on whether the time argument lies on the upper or lower branch of the contour a correlation 
function is assigned a -|- or - sign (per time argument). Below we follow a different path, which is based 
on mapping correlation function defined on the Schwinger contour on a set of alternative functions 
which obey tractable transport equations. 

Starting from the equation of motion for the field iplx) and its conjugate one can establish a hierarchy 
of coupled equations of motions for correlations functions involving increasing number of fields (in 
equilibrium this hierarchy is know as the Martin-Schwinger hierarchy [21] )• For the single particle 
propagator the equation of motion is 

Go(l)"^G(l, 1') = 5(1 ! d2dUA V{12- 34) G2(34, 1'2+), (3) 

Jc 

where Gq{1) = (i9f + V^/2m — /i)^^ and G2 are the free single particle and the two-particle propagators, 
IJ, is the chemical potential, the notation 1+ = {xi,ti + 0), the time integration goes over the contour 
and the counter-ordered delta function is defined as S{t — t') = 6{t — t') if t,t' G C+, S{t — t') = —6{t — t') 
if t, t' G C_ and 5{t — t') = {] otherwise; here C+/_ refer to the upper and lower branches of the contour 
in Fig. 121 To solve Eq. Q we need an equation of motion for the two-particle propagator which in turn 
depends on the three-particle propagator and so on. The hierarchy is (formally) decoupled by defining 
the contour self-energy as 

S(l, 3)G(3, 1') = -i ! d2dA V{12] 34) ^2(34, 1'2). (4) 
Jc 

This leads to a closed equation for the single particle propagator, which upon subtracting its complex 
conjugate takes the form 

[GS(l') - Go(l)] G(l, 1') = f d2 [G(l, 2) S(2, 1') - S(l, 2) G(2, 1')] • (5) 

Jc 

If the time arguments of the contour ordered propagators are constrained to the upper/lower branches 
of the contour we obtain the causal/acausal propagators of the ordinary propagator time perturbation 
theory 

G\l, V) = -2(r^(l)^t(i/))^ G^{1, 1') = -z(^^(l)^t(i'))^ (6) 

where T and A are the time ordering and anti-ordering operators. The fundamental difference to the 
ordinary theory is the appearance of the propagators with fixed time arguments (which can be located 
on either branch of the contour) 

G<(1, 1') = -z(^(l)^t(i')), G>(1, 1') = z(^t(i')^(i)). (7) 
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The propagators © and ((7j) are not independent 



G^(l,l') = e{U-t[)G>{l,l') + e{t[-U)G<{lA'), (8) 
= 9it[-ti)G>il,l') + e{h-t[)G<il,l'), (9) 

where 6{t) is the Heaviside step function. The equihbrium properties of the system are most easily 
described by the retarded and advanced propagators, which obey integral equations in equilibrium. 
These are defined as 

= 9{h-t[)[G>{l,l')-G<{l,l')], (10) 
G^(l,l') = 9{t[-t^)[G<{l,l')-G>{l,l')]. (11) 

There are six different self-energies associated with each propagator in Eqs. (f7j)- (fTT|) . The components 
of any two-point function defined on a time contour (in particular the single-time Green's functions and 
self-energies) obey the following relations 

:^^(i,i') -^^(1,1') = ^>(i,i') -^<(i,i'), (12) 

^^/^(l,l') = Re^(l,l')±«Im^(l,l'), (13) 

from which, in particular, we obtain a useful relation 2i Im JF(1, 1') = JF>(1, 1') — JF<(1, 1'). In practice, 
systems out of equilibrium are described by the time evolution of the distribution function which, as we 
shall see, is related to the propagator G^{1, 1'). The propagator G^(l, 1') is related to the distribution 
function of holes. To obtain an equation of motion for these functions from the equation of motion of 
the path ordered Green's function (0) we shall use algebraic relations for a convolution of path-ordered 
functions, 

C(l,l')= / rf2A(l, 2)5(2,1'), (14) 
Jc 

known as the Langreth-Wilkins rules |2S1- These rules are stated as 

C>'<(1,1') = y"c/2 [A^(1,2)5>'<(2,1') +A>'<(1,2)S^(2,1')] , (15) 

C^/^(l,l') = y"c/2 [A^/^(l, 2)5^/^(2, 1')] . (16) 

Upon applying the rule (jl5|) to Eq. (0) and using the relations (I12p and (jl3p one obtains the Kadanoff- 
Baym transport equation |5j 

[G,\l') - Re S(l, 1'), G<(1', 1)] - [ReG(l, 1'), S(l', 1)] 

= i {G>(1, 1'), S<(1', 1)} - 1 {S>(1, 1'), G<(1', 1)} . (17) 

The first term on the 1. h. side of Eq. p7|) is the counterpart of the drift term of the Boltzmann equation; 
the second term does not have an analog in the Boltzmann equation and vanishes in the limit where 
the particles are treated on the mass-shell. The r. h. side of Eq. (jl7|) is the counterpart of the collision 
integral in the Boltzmann equation, whereby S^'^(l, 1') are the collision rates. An important property 
of the collision term is its symmetry with respect to the exchange > ^ <, which means that the 
collision term is invariant under the exchange of particle and holes. Before turning to the evaluation 
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of the self-energies we briefly outline the reduction of Eq. (fTTj) to the Boltzmann's quantum kinetic 
equation. 

If the characteristic inter-collision distances are much greater than the inverse momenta of parti- 
cles and the relaxation times are much larger than the inverse particle frequencies, quasiclassical ap- 
proximations is valid. This means that the dynamics of slowly varying center-of-mass four-coordinate 
X = {xi + X2)l1 separates from the the dynamics of rapidly varying relative coordinate ^ = x\ — X2- 
One performs a Fourier transform with respect to the relative coordinates and expands the two-point 
functions with respect to (small) gradients of the slowly varying center-of-mass coordinates |2S1- Upon 
keeping the first order gradients one obtains 

z{ReG-^(p,x),G<(p,x)}p.B, +i{S<(p,x),ReG(p,x)}p.ij. = S<(p, x)G>(p, x) - S>(p, x)G<(p, x), (18) 

where P. B. stands for the Poisson bracket defined as 

{/, 9\p.B. = d^f dtg - dtf d^g - dpf drQ + d^f dpQ. (19) 

Instead of working with the functions G^''^{p,x) we introduce two new functions a{p,x) and f{p,x) 
defined by the relations, known as the Kadanoff-Baym (KB) ansatz, 

- iG<{p,x) = a{p,x)f{p,x), iG>{p,x) = a{p,x)[l - f{p,x)] (20) 

The KB ansatz is motivated by the Kubo-Martin-Schwinger (KMS) boundary condition on the Green's 
functions G^{p) = —exp[/3{uj — ^)]G^{p), where (3 is the inverse temperature, which is valid in equilib- 
rium. The KMS boundary condition is consistent with Eqs. (pUj) if we define a{p) = i [G^{p) — G^{p)] 
and identify the function f{p) with the Fermi-Dirac distribution function /i? = {1 + exp[/3(a; — /i)]}~^ 
Thus, the KB ansatz extrapolates the exact equilibrium relations (PUI) to non-equilibrium case, whereby 
the Wigner function f{p,x) plays the role of non-equilibrium distribution function which should be 
determined from the solution of appropriate kinetic equation. Eq. (jl3|) implies that in equilibrium 

^ . [G«(p) - G^(p)] ^ (21) 

which is just the ordinary spectral function, where e{p) is the free single particle spectrum. In non- 
equilibrium the spectral function need not have the form ()21|). Furthermore, the self-energies S(p) 
are functionals of the Green's functions G-^''^{p,x) and a complete solution of the problem requires 
simultaneous treatment of functions f{p,x) and a{p,x). The spectral function is determined by the 
following (integral) Dyson equation 

G^M(p)= [c.-e(p)-S«/^(p)]-\ (22) 

The level of sophistication of the kinetic equation depends on the the spectral function of the system, i.e. 
the form of the excitation spectrum. The spectral function of nuclear systems could be rather complex 
especially in the presence of bound states; for many practical purposes the quasiparticle approximation 
supplemented by small damping corrections is accurate. In this approximation j^Zl |2H] 

a{p) = 27Tz{p,x)6{uj - e{p)) --f{p,x)d^ ( , (23) 

where V stands for principal value. The first term corresponds to the quasiparticle approximation, while 
the second terms is the next-to- leading order expansion with respect to small Im S-^(p, x) or equivalently 
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small damping 7(p, x) = i x) — x)]. The wave function renormalization, within the same 

approximation, is defined as 



Note that the approximation ()23|) fulfils the spectral sum rule 



(2vr; 



A{v,x) = \. (25) 



We shall use below the small damping approximation (j25|) to establish the second and third virial cor- 
rections to the equation of state of Fermi-systems. The kinetic equation is obtained upon decomposing 
the Green's functions in the leading and next-to-leading terms in 7(p, x) ^2^ 



G>'<(p,x) =G>r(p,x)+G>j<(p,x). (26) 
Substituting this decomposition in Eq. (fTHj) one obtains two kinetic equations [211 EOl IHI] 

z{ReG'-\G<(p,x)-5S(p,x)G<(p,x)}^^ = S>(p, x)G<(p, x) - G>(p, x)S<(p, x), (27) 
i jReG^"'^, (^^^(p, x) + (5S(]9, x)G'J^](p, x)| = — ^{S^d}, x),ReG'(p, x)}p^ . (28) 

We see that the second term in Eq. (fTS)) drops out of the quasiparticle kinetic equation (j2ZI)- The 
frequency dependence of the drift term of the kinetic equation (j27j) is now constrained to have a single 
value corresponding to the quasiparticle energy; an integration over the frequency gives 

^ [S>(p, x)G'<(p, x) - G'>(p, x)S<(p, x)] . (29) 

The drift term on the 1. h. side has the familiar form of the quasiparticle Boltzmann equation; the r. h. 
side is an expression for the gain and loss terms of the collision integrals in terms of the self-energies. 
The conservation laws for particle number, momentum and energy now can be recovered from the 
kinetic equation (j29j) : e.g. integrating over the momentum we obtain the particle number conservation 
as 

The collision integrals must vanish in equilibrium. This constrains the form of the self-energies S'*'^(p, x) 
to be symmetric under the exchange > ^ <. A fundamental requirement that follows from 
the conservation laws is that the self-energies must be symmetric with respect to the interchange of 
particles to holes. In other words, the kinetic theory implies that any many-body approximation to the 
self-energies needs to be particle-hole symmetric. 



2.2 The ladder T-matrix theory 

The nuclear interactions, which are fitted to the experimental phase shifts and the binding energy of 
the deuteron, are characterized by a repulsive core which precludes perturbation theory with respect 
to the bare interaction. The existence of low-energy bound state in the isospin singlet and spin triplet 
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Figure 3: Coupled equations for the T-matrix (upper line) and the self-energy (lower line). 
The T-matrix is represented by the square, the bare interaction \^ by a vertical dashed line, 
the solid lines correspond to single-particle Green's functions. 



^5*1 — ^ Di state - the deuteron - implies further that the low-energy nuclear interactions are non- 
perturbative. The T-matrix (or ladder) approximation, which sums successively the ladder diagrams of 
perturbation theory to all orders, provides a good starting point for treating the repulsive component 
of nuclear interaction. The obvious reason is that the free-space interactions are fitted to reproduce the 
experimental phase-shifts below the laboratory energies 350 MeV and the deuteron binding energy by 
adjusting the on-shell free-space T-matrix. The contour-order counterpart of the free space T-matrix 
reads (Fig. El first hne) 

T(12; 34) = V(12; 34) + i / d5d6 V{12; 34) G(35) G(46) T(56; 34). (31) 

Jc 

The time dependence of the T-matrix is constrained by the fact that the interaction is time-local 
V{12; 34) = V{xi,X2; x^, X4)5(ti— t2)<^(i3— ^4); therefore we can write T(12; 34) = T{xi,X2, ti, X3, X4, ^3). 
For the same reason, the time-structure of the propagator product in the kernel of Eq. pip is that of a 
single two-particle propagator G{12) G(23) = G2(12;34). The retarded/advanced components of the 
T-matrix are obtained by applying the Langreth-Wilkins rule (|TH|l to Eq. (jH^. The Fourier transform 
of the resulting equation is 

T^/^(p, p'; P) = V{p, P') + P")G^'^{P"-^ P)T^'\p\ p'; P), (32) 

with the two-particle propagator 

where p± = P/2±p and the four- vector P = {P,Q) is the center of mass four-momentum. The 
remaining components of the T-matrix are given by the relations 

T>^<ip,p'■,P) = ^ J ^^T-(P,^;P)G>'<(^2)G>'<(P3)T-(^,P';P), (34) 

which can be interpreted as a variant of the optical theorem. Indeed due to the property = [T^]* 
the product T^T^ = |Tp on the r. h. side of Eq. (jM)) (we use operator notations for simplicity). At 
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the same time — T*^ = 2ilmT^ [see Eqs. ()12|) and p3|) ] which imphes that ImT^ oc G'2|7'p, where 
G2 is defined by Eq. fl^ . Thus computing T> — T< and comparing the result to Eq. (jH^ we arrive at 
another form of the optical theorem 

(2z)-iT<(p,p';P) = /M/(.;')[l-/M-/(^')]-'lmT(p,p';P) 

= gBiuj + uj')lmT{p,p;P), (35) 

where the second relation follows in equilibrium limit with gsi^) = [1 ~ exp(/5c(j)]~^ being the Bose 
distribution function. The contour ordered self-energy in the T-matrix approximation is defined as 
(Fig. El second line) 

S(l,2) = i / rf3rf4 T(12;34)G(43+). (36) 
Jc 

Since the time dependence of the T-matrix is constrained by the time-locality of the interaction, we can 
immediately write-down the two components 

^S>'<(p) = /7^T>'<f^,^;p + y') G<>{p'). (37) 



(2 



TT 



A 



where the index A stands for the anti-symmetrization of final states. Explicit expressions for the 
retarded and advanced components of the self-energy can be obtained, e.g., from the relation — = 
2zImE^ and the Kramers-Kronig relation between the real and imaginary parts of the self-energies [7j. 
Alternatively we can use the Langreth-Wilkins rules to obtain (in operator form) 

^^R,A ^ tR,Aq< + T<G^'^. (38) 

When the self-energies fl37|) are substituted in the kinetic equation ()29|) one finds the Boltzmann trans- 
port equation where the collision integrals are evaluated in the T-matrix approximation [3 1201 1221 • The 
on shell scattering T-matrix can be directly expressed through the differential scattering cross-section 

where m* is the effective mass of the particle. Thus, in the dilute limit and at not too high energies 
the collision integrals can be evaluated in a model independent way in terms of experimental elastic 
scattering cross-sections. In dense, correlated systems one needs to take into account the modifications 
of the scattering by the environment, in this case the drift and collision terms are coupled through the 
self-energies. 

Consider now the equilibrium limit. In this limit the fermionic distribution function reduces to the 
Fermi-Dirac form. The number of unknown correlations functions is reduced from two to one because 
one of the equations (PU)) is redundant. For a complete description of the system the coupled equations 
for the T-matrix and self-energy need to be solved. These are given by Eq. (j32j) where the retarded 
two-particle Green's function is now defined as 

where Q2{p+,P-) = [1 — fF{p+) ~ fpip-)] is the Pauli-blocking function. The retarded self-energy is 
given by the equilibrium limit of Eq. (jnH|) which, upon using the optical theorem ()35|). becomes 



(2vr) 



T {^—,^—;p + p)a{p)f{uj) 

R,P-P' P-P'. , /N /" du;a{p',iu) 



+2g{u + u')lmT^{P-^,P-^;p + p') J 



2n uo' — uj 



(41) 
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Eqs. flH^ . (jinj) and (PT|) form a closed set of coupled integral equations. If the interaction between the 
fermions is known these equations can be solved numerically by iteration. In the context of nuclear 
physics this scheme is known also as the Self- Consistent Green's Functions (SCGF) method 012211221 
EllESlEaEIlEHlEnilinillllliailSlliliaiini- Once the single particle Green's function (or equivalently 
the self-energy) is determined, the free energy of the system can be computed from the thermodynamic 
relation 

F = E-(3-^S, (42) 

where the internal energy is 

E = 9j ^\ + e{p)] a(p)/H^), (43) 

and the entropy S is given by the combinatorial formula 

S = 9j -^Ap) {fF{u:)\nfF{u) + [1 - fp{u)]Hl - fp{u)]] . (44) 

Here g is the spin-isospin degeneracy factor; 9 = 2 for (unpolarized) neutron matter and (7 = 4 for isospin 
symmetric nuclear matter. An important feature of the T-matrix theory is that it preserves the particle 
hole-symmetry which, as we have seen, is fundamental for the conservation laws to hold. These can be 
verified by integrating Eq. (j29j) with appropriate weights to recover the flow equations for the energy 
and momentum. Another attractive feature of this theory is that its low-density (high-temperature) 
limit is the free-space scattering theory. The latter can be constrained by experiments. The structure 
of the theory and the numerical effort needed for its solution is simplified in this limit, since instead 
of working with the full spectral function (PT|) one can approximate it with the 7(p) = limit, i.e. a 
5-function. Another interesting limit is that of low temperatures. If the damping is dropped, but the 
renormalization of the on-shell self-energies is retained (i.e. the real part of the self-energy is expanded 
with respect to small deviations from the Fermi-momentum pp) the spectral function reduces to 

a{p) = 27f Z{p)6{u - ^P)), = Pf{p - PF)/m* - /i*, (45) 

where /i* = — e(pF) + /i — ReS(pi?) is the effective chemical potential; the effective mass and the wave 
function renormalization are defined as 

* / \ -1 

m I m 



m 



(^l + ^9pReS(p)|p=,,^ , Z{p) = {l-d^ReJ:{u;,pF)l=^r\ (46) 



With these approximations one recovers the elementary excitations of the Landau Fermi-liquid theory - 
the dressed quasiparticles. Retaining the quasiparticle damping, i.e. using the small 7(p) approximation, 
Eq. pHjl . leads to virial corrections to the quasiparticle pictures. We shall discuss these corrections in 
more detail in Subsec. 12.2.41 

The T-matrix approximation to self-energy leads to a model which satisfies the conservation laws 
(it is said that the model is conserving). In addition to being conserving any model that is based on an 
certain approximation to the self-energy needs to be thermodynamically consistent. The thermodynamic 
consistency refers to the fact that thermodynamic quantities like free energy or pressure computed from 
different expressions agree. An example is the Hugenholtz-van Hove theorem j32j, which relates the 
single particle energy at the Fermi-surface to the binding energy Eb at the zero temperature 

e{PF) = - + Eb, (47) 
P 
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Figure 4: Left panel. Dependence of the experimental scattering phase shifts in the ^5*1, 
^P2, ^-02, and ^Di partial waves on the laboratory energy. Right panel. The dependence 
of the critical temperatures of superfluid phase transitions in the attractive channels on 
the chemical potential. The corresponding densities are indicated by arrows. The critical 
temperatures Tc are computed from the T-matrix instability [1^ . 



where the pressure is defined as P = p^dpEs- Another example is the equivalence of the thermodynamic 
pressure defined above and the virial pressure, the latter being the pressure calculated from the energy 
momentum tensor. 



2.2.1 Pairing instability and precursor phenomena 

We have seen that in the low-density and high-temperature domain the T-matrix theory is well defined 
in terms of free-space parameters and it can be used at arbitrary temperatures and densities. In the 
opposite limit of high densities and low temperatures its validity domain is restricted to the temperatures 
above the critical temperature Tc of superfluid phase transition. The physical reason is that at Tc 
there appears a bound state in the particle-particle channel - the Cooper pair. This has far reaching 
consequences, since the onset of macroscopic coherence implies that the average value of correlation 
function {il:{x)ip{x)) ^ 0, which requires a doubling of the number of Green's function needed to 
describe the superfluid state. At temperatures T > Tc the T-matrix is strongly enhanced for particle 
scattering with equal and opposite momenta and it diverges at T = Tc. 

Partial wave analysis of the nucleon-nucleon scattering allows us to identify the attractive channels 
(which feature positive phase shifts). The critical temperature in each channel is determined from the 
condition that the T-matrix, Eq. (jH^ . develops a pole for parameter values P = {fl, P) = (2/i, 0) [7| 
I2H1 EHl Eni Eni • To illustrate this feature assume a rank-one separable interaction V{p,p') = x{p)x{p') 
and the quasiparticle approximation. The solution of Eq. ()32j) . which parametrically depends on the 



14 




Figure 5: The cross section in neutron matter as a function of relative energy of particles 
[E = £^lab/2) for temperatures T = 10 MeV {left panel) and T = 4 MeV {right panel) and 
densities p = (solid line), p = po/5 (dashed line), p = po/2 (dashed-dotted line) and p = po 
dotted line. 



chemical potential and the temperature is 

(48) 

At the critical temperature both the real and imaginary parts of the expression in braces vanish; the zero 
of the real part determines the critical temperature Tc. Fig.|3]shows the neutron-proton scattering phase 
shifts which are relevant for the pairing pattern in the isospin symmetric nuclear matter (left panel) 
and the associated critical temperatures (right panel) determined from the T- matrix instability |49j . 
For isospin symmetric systems the most attractive channels are the tensor channel ^Si Di and 
the channel where only the neutrons and protons interact. For small isospin asymmetries, which 
correspond to a = {pn — Pp)/{pn + Pp) < Q^c — 0.1, where p„ and pp are the neutron and proton densities, 
the mismatch in the Fermi-surfaces of neutrons and protons suppresses the pairing [301 1^ 13^ EH 
l55| IHEl inZl IHHl inn] • For large asymmetries typical for compact stars the pairing is among same isospin 
particles in the ^5*0 and channels. Because of the smallness of the charge symmetry breaking effects, 
the critical temperatures in the ^5*0 and shown in Fig. |3] channels are representative for neutron 
star matter as well (note however that the relation between the density and the chemical potential 
changes) jSOllS]- Some models of neutron star matter which predict kaon condensation at high densities 
feature isospin symmetric nucleonic matter, in which case the high-density D-wave neutron-proton 
paring will dominate the P-wave neutron-neutron pairing [HU]. The scattering characteristics of the 
system such as the phase-shifts and the scattering cross-sections are affected by the pairing instability, 
since these are directly related to the on-shell T-matrix. The phase-shifts in a pairing channel change 
by 7r/2 at the critical temperature when the energy equals 2p. According to the Levinson theorem. 



T{p,p\P) = V{p,p') 



1 - 



(2^ 



x\p)G2{p,P) 



15 



this corresponds to the appearance of bound state (Cooper pair). For many practical apphcations the 
cross-section is the relevant quantity. According to Eq. (|39|) the cross-section being proportional to 
the T-matrix will diverge at j2Hl UHl in21 IHSl IHH IHHl inHj ■ The precursor effect of the superfluid phase 
transition on the neutron-neutron scattering cross-section in the low temperature neutron matter is 
shown in Fig El ^3 • The cross-section develops a spike for lower temperature as a precursor of the 
onset of superfluid in neutron matter in the ^5*0 interaction channel. The largest enhancement is seen 
for the density p = 0.5po which is closest to the maximum of the critical temperature as a function 
of density. The above precritical behavior of the cross-section has a significant effect on the transport 
and radiation processes in matter, for example, it could lead to a critical opalescence in the transport 
phenomena. 



2.2.2 T-matrix theory in the superfluid phase 

We have seen in the previous section that at the critical temperature of superfluid phase transition the 
two-body scattering T-matrix develops a singularity, which is related to the instability of the normal 
state with respect to formation of Cooper pairs; this is manifested in pole of the two-body T-matrix 
when the relative energy of interacting fermions is twice their chemical potential. Thus, the T-matrix 
theory described above breaks down at the temperature T = T^. A T-matrix theory appropriate for 
temperatures below T^ can be formulated in terms of the normal and anomalous Green's functions [7j. 
To account for pair correlation we represent each Green's function in the Keldysh-Schwinger formalism 
as a 2 X 2 matrix in the Gor'kov space: 

ir(. = ( F^p{x,x') \ / -i{T^^{x)^l{x')) {^^{x)^p{x')) \ 

where Gap{x^x') and i^^^(a;, x') are referred to as the normal and anomalous propagators. The 4x4 
matrix Green's function satisfies the familiar Dyson equation 



C„^(a;,x') = Gl%{x,x') + fd^x''d*x'''G!%{x,x''')^^s{x''\x'')G!sf,{x'',x') 

7,<5 



(50) 



where the free propagators ^°^(x,x') are diagonal in the Gor'kov space. We consider below uniform 
fermionic systems; the propagators now depend only on the difference of their arguments due to trans- 
lational symmetry. A Fourier transformation of Eq. (j^U)) with respect to the difference of the space 
arguments of the two-point correlation functions leads to on- and off-diagonal Dyson equations 

Ga^ip) = Goa/3(p) + Goa^/ip) [^^s{p)Gsp{p) + A^s{p) Fdsp{p)] , (51) 
FKpip) = Goa^{-p)[A\,{p)Gsp{p) + i:,s{-p)FAsf3{p)], (52) 

where p is the four-momentum, Goapip) is the free normal propagator, and Tia/B^p) and Aapi^p) are 
the normal and anomalous self-energies. Summation over repeated indices is understood. Specifying 
the self-energies in terms of the propagators closes the set of equations consisting of (j5ip and (j52j) and 
their time-reversed counterparts [Ga/sip) and i^a^(p)]. The particle-particle scattering in the superfluid 
state is described by three topologically different vertices shown in Fig. IHl -nZl- We write out the explicit 
expression for the retarded components of the T-matrix in operator form [their form in the momentum 



16 



X = xx*)o<*xx 

Figure 6: The ladder series in a superfluid fermionic system. The dashed hnes correspond to 
the driving interaction, single arrow lines to normal propagators and double arrow lines to 
anomalous propagators. The filled square is the counterpart of the unpaired state T-matrix, 
the shaded and empty squares are specific to superfiuid state. 



space in identical to Eq. 

T(i) = V[1 + SggT^^'> + SfgT'^^^ + SffT'-^^] , (53) 
T(2) = V [SffT^^^ + SggT^^^ + SfgT^^^] , (54) 
T(3) = V [SfgT'-^^ + SggT^^'' + SffT^^^ , (55) 

where the two-particle retarded propagators are defined as 

Sgg = G>G> -G<G<, Sgf = G>F> -F<G<, Sgf = F>F> - F<F<, (56) 

and their structure in the momentum space is given by Eq. (j33|) . To close the system of equations we 
define the normal and anomalous self-energies for the off-diagonal elements of the Schwinger-Keldysh 
structure 

,S>'< = t(i)>.<G'<'> + t(2)>'<F<'>, (57) 
zA>'< = t(=^)>.<f<.> + t(2)>.<G'<'>, (58) 

which are shown in Fig. [7| The retarded components of the self-energies which solve Dyson equations 
dSU and can be constructed from Eqs. (jFTj) and via the dispersion relations, e.g., 

^ r ^ A^M^AIM. (59, 

with an analogous relation for S^(co'). The system of Eqs. (|^ - (f^ can be used to derive the collective 
excitations of the superfluid nuclear matter in particle-particle channel. In the mean field approximation 
the self-energies decouple from the T-matrix equations (j53p -(|55 p and the secular equation determining 
the frequencies of collective modes for vanishing center-of-mass momentum is 

[A + (4A2 - LJ^)B] [A - {AA^A + uj^)B + 2A^uj^B^] = 0, (60) 
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Figure 7: The normal (filled circle) and anomalous (empty circle) self-energies in the su- 
perfluid T-matrix theory. The self-energies couple through the T matrices appearing in 
Fig. El 



where, assuming that the pairing interaction can be approximated by a constant and the integrals 
regularized by a cut-off, one finds 

B{u;) = X r ^-^MM , /! . , A = l + X r^J^[l-2Me)], (61) 

where A is an effective coupling constant, A is an ultraviolet cut-off. The second equation is the 
stationary gap equation, i.e. A = 0; the secular equation then leads to the solutions Bioj) = 0, 
B[lj) = (2A^)~^, uj"^ = and uj = ±2 A. Among the first two non-trivial conditions the second one does 
not have a solution for weakly coupled systems A <^ 1 and the collective modes are determined by the 
secular equation B{u!) = 0, whereby the real part of the solution determines the eigenmodes and the 
imaginary part their damping. 

2.2.3 Three-body T-matrix and bound states 

Up to now we were concerned with the correlations described by the two-body T-matrix. The properties 
of dilute fermions or cold Fermi-liquids (the latter are characterized by a filled Fermi-sea) are well de- 
scribed in terms of two-body correlations between particles or quasi-particle excitations. However, the 
three-body correlations, which are next in the hierarchy, are important under certain circumstances. We 
turn now to the three-body problem in Fermi-systems within the formalism developed in the previous 
sections. As is well known, the non-relativistic three-body problem admits exact free space solutions 
both for contact and finite range potentials [HHl 111] ■ Skorniakov-Ter-Martirosian-Faddeev equations 
sum-up the perturbation series to all orders with a driving term corresponding to the two-body scatter- 
ing T-matrix embedded in the Hilbert space of three-body states. The counterparts of these equations 
in the many-body theory were first formulated by Bethe [TO] to access the three-hole-line contributions 
to the nucleon self-energy and the binding of nuclear matter (Bethe's approach is discussed in Subsec- 
tion EISI)- More recently, alternative forms of the three-body equations in a background medium have 
been developed that use either an alternative driving force (the particle-hole interaction or scattering 
T-matrix) jSHl IZH IZ21 ES] or/and adopt an alternative version of the free-space three-body equations, 
known as the Alt-Grassberger-Sandhas form [Z2 ESI EE] ■ 

The resummation series for three-body scattering amplitudes can be written down in terms of the 
three-body interaction V as [22] 

T = v + vgv = v + vgor, (62) 
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Figure 8: Coupled integral equations for the three-body T-matrix (first line) and the single 
particle propagator (second line). The shaded vertex stands for the amplitude T^*-*, the 
empty vertex stands for a channel Tkj matrix. The second line shows the Dyson equation for 
the single particle Green's function including the contribution from three-body scattering. 
The dots stand for the contribution from the two-body scattering shown in Fig El 



where Qo and Q are the free and full three-particle Green's functions (we use the operator form for 
notational simplicity; each operator, as in the two-particle case, is ordered on the contour). If the 
three-body forces that act simultaneously between the three-particle are neglected, the interaction in 
Eq. (j62|l is simply the sum of pairwise interactions: V = V12 + V23 + V13, where Va/j is the interaction 
potential between particles a and /?. The kernel of Eq. is not square integrable: the potentials Vap 
introduce delta-functions due to momentum conservation for the spectator non-interacting particle and 
the iteration series contain singular terms (e.g., of type Va/sGo'^afs to the lowest order in the interaction). 
The problem is resolved by summing up the ladder series in a particular channel (specified by the indices 
a, 13) to all orders This summation defines the channel T^/j-matrix, which is essentially the two- 
body T-matrix embedded in the Hilbert space of three-particles states 

= Vq,/3 -|- Va/3 Go %yf3- (63) 

The three-body T-matrix can be decomposed as T = T^^^ + T^"^^ + T*^^) , where 

T(") = Vp^^ + Vp-,GoT, (64) 

and a(3'-) = 123, 231, 312. Now, Eqs. ()63p and ()64|) are combined to eliminate the interaction terms 
Va/3 and one is left with three coupled integral equations for T^") (a = 1, 2, 3) 

T(-) = T^, + Tf,, g, (T(^) + T(^)) , (65) 

where the driving terms are the channel T-matrices. The new equations are non-singular Fredholm 
type-II integral equations. Note that their formal structure is identical to the Faddeev equations in the 
vacuum [Ml; however their physical meaning is different. To see the physical content of Eqs. (jU^ we 
need to convert the contour ordered equations into equations for the components (so that the KB ansatz 
(|2Up can be applied) and to transform them from the operator form into momentum representation. 
Proceeding as in Subsection 12.11 the retarded component of Eq. reads 

r^(i)(t, t') = T^f (t, t') + l[ T^^'\t, t) + r^^'\t, t)] g^it, t")%^{t", t')dtdt", (66) 
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were we used the time-locality of the interaction and omitted the momentum arguments of the functions 
(for the explicit expressions see ref. [29j). Next, to apply the KB ansatz we need to specify the particle 
hole content of the three-body T-matrix, i. e. assign each incoming/outgoing state a particle or a hole. 
Figure |H1 shows the Feynman diagram for the three-body T-matrix where all the incoming (outgoing) 
states are particles (holes). The remaining three-body T-matrices are obtained by reverting the direction 
of the arrows in the diagram. Depending on the particle-hole content of the three-body T-matrix (in 
the sense above) the intermediate state retarded Green's function is 

G>G>G>{h,t2)-i>^<) (3j9) 

G<G<G<{hM)-{>^<) i3h) 

where p and h refer to particle and hole states, the brackets [e. g. (2ph)] indicate the particle-hole 
content of Green's function; for simplicity the time argument in a product is shown only once. The 
short hand > ^ < stands for a term where all the G^ and G^ functions are interchanged. Upon applying 
the KB ansatz and Fourier transforming Eq. (jMj) one finds 



(68) 



where the the four- momentum space is spanned in terms of Jacobi coordinates, K = Pa+Pis+P-y, ka/s = 
{pa — P/3)/2, q-y = {pa + Vp)/'^ ~ 2^7/3, the center-of-mass energy VL = Kq and 

QsiPmPlS^Pj) = a{pa)a{p(i)a{pj) {[1 - fpiPaW - friP/sW - fpiPj)] - fF{Pa)fF{P(})fF{Pj)} ■ (69) 

This form of three-body equation incorporates off-mass-shell propagation if the spectral function is 
taken in the form (PT|) . The quasiparticle (on-mass-shell propagation) limit follows by using in Eq. 
(jnH). Thus, the many-body environment modifies the three-body equation in a twofold way: first, the 
single-particle spectrum is renormalized in the resolvent of Eq. ()68|1 [which becomes explicit after taking 
the quasiparticle limit], second, the intermediate state propagation is statistically occupied according 
to Eq. (jnH). The limit Q — >■ 1 and a (a;) = 27i6{u! — e(p)) recovers the original Faddeev equations. The 
single-particle self-energy obtains contributions from the three-body T-matrix, which is shown in Fig. |H| 
In the case of (3p) scattering T-matrix the three-body self-energy is written as 

©>'<(pi) = / \ri]l:ik2s,quK) + ri%ikrs,q2,K) + Ti%ik,,,qs,K)]^^ (70) 



where we use the notation j = d'^p/ (27r)^. The optical theorem relates the T^'^ components of the 
three-body T-matrix to the retarded component given by Eq. ()68|1 : 

T<{k2s,quK3^q[;K) = / T^'^^ ""{k^s, q^ k,,, q^, K)G<ip,) G<{p,) G<{p,) 

xT(«^(fc45, 96, K„ q[; K){2n)^ 5\K - p, - p, - p,). (71) 

To illustrate the usefulness of the three-body equations discussed above we show in Figure IHl the binding 
energy of a three-body bound state (triton) in nuclear matter as a function of the inverse temperature 
for several values of density, n, measured in units of a~'^, where a = 5.4 fm is the neutron-proton 
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Figure 9: Dependence of the three-body and two-body bound state energies on the inverse 
temperature for fixed parameter values na^, where n is the density and a the neutron-proton 
triplet scattering length. 



triplet scattering length. The asymptotic free-space value of the binding energy in this model is E^b = 
7.53 MeV. The binding energy was obtained from the solution of the homogeneous counterpart of 
Eq. (j68|l in quasiparticle hmit assuming free single particle spectrum jUj. In addition we show the 
temperature dependence of the deuteron energy E2b{P) obtained within the same approximations from 
the homogeneous counterpart of Eq. (jHSI)- The continuum for the break-up process 3B 2B + A^, 
where N refers to nucleon, is temperature/density dependent as well and is found from the condition 



2.2.4 The quantum virial equation of state 

The equation of the state of a Fermi-system characterized by small damping (long-lived, but finite 
life-time quasiparticles) can be written in the form of a virial expansion for density j221 |2H1 EHj 

^ [aiE,p)ME) + b{E,p)gB{E) + c{E,p)ME)] , (72) 

where the virial coefficients a{p), b{p), and c{p) are the one-, two- and three-particle spectral functions. 
Below we show that the virial coefficients b{E,p) and c{E,p) can be written entirely in terms of the 
two- and three-body T-matrices and their derivatives. In the dilute limit the on-shell T-matrices are 
related to the scattering phase-shifts; since the damping in this limit is small a direct relation between 
the scattering observables in free-space and the equation of state can be established. We have seen that 
in the small damping limit the spectral function can be approximated by Eqs. (j^Hj) and ()24|) to leading 
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order. Thus our starting point is the expression for the density of the system which we write as 

/d'^p f d^p 

+ j ^^«M [/f(^) - /^(e)] • (73) 

The first term is the contribution from the "uncorrelated" quasiparticles (note that the notion of 
quasiparticle already requires correlations which renormalize the single particle spectrum, however the 
quasiparticles are still characterized by a sharp relation between the energy and the momentum as the 
ordinary particles). The second term is the correlated density which, upon using Eqs. and (j^ . 
becomes 

ncorr(/5, ^^) = - ! -^Mp)d. [/f(^) - fF{e,)] . (74) 



Let us first evaluate this expression neglecting the three-body correlations. The damping can be written 
as 7(p) = z[S^(j9) — S^(p)] = — 2ImS where the self-energies are given by Eqs. (|n7|). Using the optical 
theorem for the two-body T-matrix we obtain 

7(Pi) = 2 j ^«(^2) [^?B(cui+c.2) + M^2)]Imr^(^^^,^^^;Pi+P2), (75) 

where the spectral function can be taken in the quasiparticle approximation at the order of interest. 
Substituting the damping in Eq. (fflj) and using the identity [gB^E) + fril^i)] [fpiE — 002) — fri^i)] = 
[gsiE) — gsi^^i + ^2)] [1 — fpi^i) — /f('^2)] , we recover the second term of the expansion (ffHjl with the 
second (quantum) virial coefficient (23 12H] 



(76) 



where R^{E) = [E — e{pi) — e{p2) + irj] ^ is the two-particle resolvent. For systems which support 
bound states in the free space the second virial coefficient obtains contributions both from the negative 
energy bound states and the continuum of scattering states. The bound states appear as simple poles 
of the two-body T-matrix on the real axis. The scattering states can be characterized by the phase- 
shifts in a given partial wave channel, after the two-body scattering T-matrix is expanded into partial 
waves. The phase shift is defined simply as the phase of the on shell complex valued matrix To,{p,p; E = 
e{p) + e{p)) = \Ta{p,p, E)\exp {6a) where a = TSJLL' specifies the partial wave channel in terms of 
total spin S, isospin T and angular L and total J momenta. The second virial coefficient can now be 
written in terms of the bound state energies Ep (where (3 = 1, 2... enumerate the poles of the T-matrix) 
and the scattering phase shifts 

b{e,„E) = 2nJ2H^ - Ep) + 2sm'6aiE)^^^, (77) 

P a 

where are channel dependent constants. In the non-degenerate limit one recovers the classical Beth- 
Uhlenbeck formula [75] . 

The third virial coefficient is obtained by including the contribution to the damping from the three- 
body processes C5 = — = — 2Im(3, where the self-energies are defined by Eq. dZOl). Now the 
correlated density is written as 

n.oAP,f^) = - I^Mp) + ^(P)]^- (^7^) f^^^""^ ~ ^^^^^^^ • ^^^^ 
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The damping (i5(p) is expressed in terms of the T^'^ matrices which in turn can be related to the 
retarded component if we use the optical theorem obeyed by the three-body matrices. The off-shell 
form of the optical theorem reads 

lmr''{k2sqi,k,,q,;K) = I r''^''\k,,q,, k,,q,; K)\G<{p,)G<{p,)G<{p,) 

-G>{p,)G>ip,)G>ipe)] r^^^\k,,q,, k^squ K)5' {K - p, - p, - p,) , (79) 

where the proper momenta pi {i = 1,2,3) and Jacobi momenta kij qk transform into each other according 
to the rules given after Eq. ()68|) . Since we seek corrections that are first order in damping, the T-matrix 
in expression (|78p can be taken in the quasiparticle approximation; the on-shell optical theorem implies 
then 

^^(^2391,^45^6;^) = 2ifFiuJi + uj2 + uj-i)lmT^{k23qi,ki5qG;K), (80) 

and an analogous expression for with the replacement fpi'^) — > 1 — fri'^)- Evaluating the three- 
body damping C5(p) with the help of on-shell optical theorem one arrives at the third quantum virial 
coefficient ^29j 

c{ep,,E) = 2 / Q3{ep„ep^,ep^) lmT^{k23qi,k23qi; K)—Ren{E) 

+ Im7^(E)^Rer(fc23gl,fe23gl;i^)], (81) 

and where IZ = [E — Sp^ — Sp^ — Sp.^ + iri]~^ is the three-particle resolvent. The third virial coefficient 
can be decomposed into scattering and bound-state contributions in analogy to the two-body case. 
Complications arise in attractive systems where apart from the three-body bound states one needs to 
take into account the break-up, recombination and rearrangement channels which are absent in the 
two-body case. The knowledge of the virial expansion (f?^ completely specifies the equation of state of 
the system; the pressure can be computed from the Gibbs equation 

p{fi,f3)= r rf/i'n(/i',/5). (82) 

J — oo 

The common form of the equation of state p{n, (3) is obtained upon eliminating the parametric de- 
pendence on the chemical potential /i. The theories based on the second virial coefficients smoothly 
interpolate between the classical gas theory at low densities and high temperatures and the Bruckner- 
Bethe-Goldstone theory at low temperatures and high densities j2H] • The effects of the third quantum 
virial coefficient on the equation of state of nuclear matter have not been studied to date. 



2.3 The Bruckner- Bethe-Goldstone theory 

In a number of cases it is more convenient to evaluate the ground state energy, or at finite temperatures 
the thermodynamic potential, directly instead of first obtaining the Green's function from the Dyson 
equations and then calculating the thermodynamic quantities. The Brueckner-Bethe-Goldstone (BBG) 
theory evaluates the ground state energy of nuclear matter in terms of certain diagrammatic expansion 
of the energy, which has two important ingredients: (i) the effective interaction is built up from the bare 
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Figure 10: The lowest order diagrams of the BBG theory for the direct {left) and the 
exchange {right) contributions to the energy. The square denotes the i^'-matrix, the sohd 
hues correspond to the hole propagators. 



nucleon-nucleon force by summing the ladder diagrams into an effective interaction; (ii) the perturbation 
expansion is organized according to the number of independent hole lines in all topologically non- 
equivalent linked diagrams [ZHl EH lEB El ESI EH ■ 

The diagrams describing the perturbation series for the energy have the form of closed loops. Only 
the connected (linked) diagrams contribute, i.e. those diagrams which have the property that by starting 
at a vertex one can return to the same vertex moving along all the interaction and propagator lines. 
The diagrammatic rules are the same as those of the ordinary Feynman perturbation theory (except of 
the overall constant factor) 0. The lowest order diagrams of the BBG theory are shown in Fig. ITUl 
The energy of nuclear matter including the contributions of the lowest order diagrams is written as 

where it is understood that at small temperatures the Fermi-functions are approximated by the step 
functions /f(p) = d{pF — p)', the superscript (2) indicates that the hole line expansion is carried out up 
to the terms of second order. The interaction is approximated by the i^-matrix (We use here the term 
/T- matrix instead of the G- matrix to avoid confusion with Green's functions). The i^'- matrix sums the 
ladder diagrams, where the driving term is the bare nucleon-nucleon interaction 

K{p, p'- p, p') = V{p, p'- p, p') + / ^V{p, p'- g, q') - ^'^f 7 ^{"5^] K{q', q- p, p') . (84) 

J {2-n:)^ uj - e{q) ~ e{q') + 17] 

In the intermediate state the i^-matrix propagates two particles; the hole- hole propagation oc f f{(i) f f{(i') 
which appeared in the T-matrix defined by Eq. (jH^ is absent here. As a consequence the i^-matrix 
does not develop a singularity at the critical temperature of superfluid phase transition and is well 
defined at zero temperature. It is clear that to obtain the equation for the i^-matrix from Eq. (j32j) the 
quasiparticle limit (j45p must be taken with the wave-function renormalization Z{p) = 1. Note that the 
effective interaction entering the BBG expansion is real, therefore the +iT] term in Eq. (j84p is commonly 
dropped and the integration is treated as principal value integration. Since the perturbation expansion 
is now carried out for a macroscopic quantity, it is not obvious what the single particle energies e{q) in 
Eq. (jH^ represent. This ambiguity leads to several choices of the single particle spectrum, one possible 
form being 
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a 



b 



c 



Figure 11: The lowest order three- hole line contributions to the energy of nuclear matter: 
(a) the third-order bubble diagram, (b) the third order ring diagram, (c) the bubble diagram 
with an auxiliary potential insertion. The particles above the Fermi sphere propagate from 
the left to the right, the holes - from the right to the left. The dots denote the iiT-matrix, 
the shaded vertex - the auxiliary potential U. 



where U{p) is called the auxiliary potential. (The term arises from the rearrangement of the original 
Hamiltonian H = T+V = T+U+6V, where 6V = V — U is small, T and V are the kinetic and potential 
energies). The so-called 'gap choice' keeps the auxiliary potential for the states below the Fermi surface, 
which leads to a gap in the spectrum at the Fermi-energy; the 'continuous choice' keeps this potential 
both for the particle and the hole states. Another definition arises upon using the Landau's Fermi-liquid 
theory, where the quasiparticle energy is defined as the functional derivative of the total energy with 
respect to the occupation e{p) = 6E/6fp{p) j|85j 



./^iA-(p..-;.„oi-/.(.)/.(.') ,/':;f;':,;-!'',,, (s^) 

The last term, known as the rearrangement term, guarantees that the quasiparticle energy is in fact 
the energy needed to extract a particle from the system. Returning to the BBG theory, it should be 
noted that the choice of the single-particle spectrum (i.e. the self-energy) specifies the set of diagrams 
that are already included in the Green's functions from which the closed diagrams for the energy are 
constructed; and it is a matter of convenience which building blocks are chosen as fundamental. The 
gap choice implies that the self-energy insertions for the particles are treated explicitly by grouping 
them into the higher order clusters. Thus, the particles and the holes are treated asymmetrically both 
in obtaining the effective i^-matrix interaction and in defining the single particle energies. Note that 
in contrast to the T-matrix theory where the self-energies S^'^(p) are defined symmetrically, the BBG 
theory breaks this symmetry. 

Given the form of the single particle spectrum, the next natural question is the organization of the 
diagrams in an expansion which has reasonable convergence properties. The BBG theory identifies such 
expansion parameter and organizes the diagrams order by order in this small parameter according to 
the number of independent hole lines (i. e. the number of hole lines that remain after the momentum 
conservation in a given diagram is taken into account). A diagram with i independent hole lines is of 
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the order of ^ where the k parameter (pair excitation probabihty) is defined as 

K = n I \^{x) ~(l){x)\^d^x, (87) 



where ip{x) is the perturbed pair- wave function satisfying the Schrodinger equation Kef) = Vip, while 
(f){x) is the uncorrelated wave function and n is the density. The integral extends to the surface where 
wave-function restores its free-space form and defines an 'interaction' volume ~ 4'7rA^/3, where A is the 
hard-core radius. Thus, we see that k is essentially the ratio of the volumes occupied by the hard-core 
interaction and a particle. At the saturation density of nuclear matter k is of the order 15%, which 
allows one to estimate the error introduced by neglecting an i-hole line diagram to the potential energy 

AEpot(z) = AEp,t{2)K'-\ 

The expansion clearly breaks down when the interparticle distance is of the order of the hard-core radius 
of the potential. 

The computation of the next-to-leading order three-hole line diagrams is complicated by the fact 
that the scattering problem of three particles in the nuclear medium needs to be solved. The appropriate 
equations are due to Bethe and are known as Bethe-Faddeev equations jZOj- These equations are the 
counterparts of the three-body Faddeev equations in the free space, which take into account the infiuence 
of the background medium. The lowest order three-hole line diagrams are shown in Fig. ^2 ,86, 87j. 
The particles above the Fermi sphere propagate from the left to the right, the holes - from the right to 
the left. The i^'-matrix is represented by a dot, since the BBG theory assumes K to be local in time. 
In the case where the spectrum is chosen to have a gap, one needs to evaluate only the diagrams a and 
b in Fig. ^2 while in the case of a continuous spectrum the diagram c in Fig. II II should be evaluated as 
well; here the shaded vertex is an insertion of the auxiliary potential U. (The rationale behind the gap 
choice is the cancellation of this type of diagrams in the BBG expansion, albeit, such a choice requires 
evaluation of the three-hole line diagrams, contrary to the continuous choice which is well converged at 
the two-hole line level) j88 | 1 ^ IHTTj. The contribution to the energy from the three- hole line terms can 
be written, using for simplicity the operator notations, as 



- ¥^ (I) 



wherex = P123 + -P132 and the three-body /C matrix satisfies the Bethe-Faddeev equation 

/C = Kax Kb{1 + x)- Kax (^) /C. (89) 



The first term in Eq. (|89|) is the sum of the direct (oc 1) and exchange (oc x) diagrams shown in Fig. 111! 
the second term corresponds to the so-called higher order diagrams. Here the permutation operator 
is defined as -PIP3P1P2) = IP1P2P3) + \P2P3P1), the index A indicates that the spectator third particle 
is above the Fermi-surface and the index B indicates that there is non. This implies that the matrix 
elements of the energy denominators in the three-body basis (here q denotes the spectator particle) are 
defined as 

_ J E{q) + E{k) + E{p) - Us (if q > qp, case A), 
^^/^ \ E{k) + E{p)-uj2 (if q<qF, case B), ^''''^ 
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Figure 12: Le/( panel: The three-hole hne contributions to the energy per particle of the 
isospin symmetrical nuclear matter as a function of the Fermi wave vector. Displayed are 
the contributions from the direct [dashed line), exchange [dashed- dotted line), higher order 
[dotted line) diagrams and their sum E^^^ /A [solid line). The heavy lines marked with squares 
show the results obtained with the Reid potential |HZ1 HH] the light lines - with the Argonne 
V14 potential jHHlEni- Right panel: The energy per particle of isospin symmetrical nuclear 
matter within the BBG theory. E^'^^ [dashed line) is the contribution which includes the 
two-hole line diagrams only; E^'^^ + E^^'' include in addition the three-hole line contributions; 
the labeling of curves and references are the same as in the left panel, however the results 
for refs. [88, 89 for the Argonne V14 potential are shown for the continuous spectrum. The 
empirical saturation point is shown with error bars. 



where cus = E[pi) + E[p2) + E[p^) and uj2 = E[pi) + E[p2), where E[p) is the single particle energy. 
The action of the three-body Pauli operator is written as 

Qs\q,kp) = [l-fp[k)][l-Mp)], fp[p)^9[\pp-p\)^ (91) 

which implies that the particles in the two-body subspace must be outside the Fermi-sphere, while the 
propagation of the third, spectator particle, is not restricted by the Pauli principle. The BBG Pauli 
operator should be compared to Eq. ()69|) which is the most general form of a three-particle Pauli 
operator for intermediate particle propagation which preserves the particle-hole symmetry. 
Fig. 1121 shows the various contributions to the three-hole line energy 

^(3) 

per particle for several 

densities parameterized in terms of the Fermi wave vector [n = 2/c|,/37r^ in symmetrical nuclear matter). 
An important feature is the mutual cancellation of the positive contribution from the direct term and 
the negative contributions from the exchange and higher order terms. The differences between the 
results of refs. jHZl EO] and [SHI IMI shown in the left panel Fig. [T21 are due to the differences in the 
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Reid and Argonne V14 potentials. The differences in the right panel are mainly due to the choice of 
the spectrum - gaped spectrum in the first case and continuous spectrum in the second. The three-hole 
line contribution to the energy leads to a saturation curve of nuclear matter, which predicts a binding 
energy that is consistent with the empirical saturation point (Fig. ^1 right panel). The minimum 
of the saturation curve lies at densities that are larger than the empirically deduced one - the missing 
ingredient is the three-body forces. The convergence of the hole-line expansion in the case of continuous 
spectrum is faster than in the case of the gaped spectrum; in the first case the two-hole line expansion 
provides a satisfactory results within the errors which are introduced by ignoring additional physics, 
such as three-body forces and relativistic dynamics (these two aspects of many body problem cannot 
be disentangled in general). 

2.4 Relativistic T-matrix theory 

In describing the nuclear phenomenology within relativistic theory two distinct approaches are possible: 
the phenomenological approach starts with a meson-baryon Lagrangian whose parameters are fitted to 
reproduce the known empirical properties. A typical set is the binding energy at saturation Eb — —16.0 
MeV, saturation density po = 0.16 fm~^, compression modulus K ~ 300 MeV, symmetry energy 
Es ~ 30 MeV (see Subsec. below) and effective nucleon mass at saturation m%- = O.Srriiy, where rriiy 
is the bare nucleon mass ^T]. The microscopic approach constructs first the free-space scattering T- 
matrix from a one-boson-exchange potential, which fits to the scattering phase-shifts and the deuteron 
binding energy; given the free-space interaction a many-body scheme is applied to describe the physics 
in matter. These models are then extrapolated to the large densities (and temperatures) to describe 
the properties of matter under stellar conditions. This bottom to top approach (with respect to energy 
scales) should be contrasted to the top to bottom approaches that attempt to constraint the form of the 
nucleon-meson Lagrangian and the couplings by the symmetries of the underlying fundamental theory 
- quantum chromodynamics (QCD). The models that incorporate the chiral symmetry - the dynamical 
symmetry of strong interactions - are based on low-momentum expansions of chiral Lagrangians; the 
usefulness of chiral models for treating dense hadronic matter, where momenta are generally not small 
compared to other relevant scales (e.g. Fermi-energies) is unclear. However, chiral models are useful in 
treating the meson nucleon interactions in matter; for example, these have been used extensively in the 
studies of the kaon nucleon interactions in matter [HI] (see Subsec. 12. 81 for a discussion and Subsec. 12.101 
for further references). 

2.4.1 Dyson- Schwinger equations and mean field 

The elementary constituents of the relativistic models of nuclear matter are the meson and the baryons, 
whose interaction can be described by a model Lagrangian 

= -9s'ip'ip(ps - igpv^T, ^(5/.0Pv) - 5'yV'7M^0v' - ^fvi'T, "^{d^, (92) 

where the gs, gv, fv and gpv are the coupling constants of the nucleon fields to the meson fields 0, 
the indices S, V, and PV refer to scalar, vector and pseudovector couplings. Table [T] lists the (non- 
strange) mesons, their quantum numbers and typical values of meson-nucleon couplings. The a meson 
is believed to represent the two-pion exchange contribution to the interaction within the one-boson- 
exchange models. Chiral symmetry of strong interactions allows the presence of self-interacting meson 
terms in Eq. ()92|) which we neglect for simplicity. The Euler-Lagrange equations for the baryon and 
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meson 


TT 




P 


a 


V 


5 


4> 


Mass [MeV] 


139 


784 


764 


571 


550 


962 


1020 


7, J'^ 


1,0- 


0,1- 


1,1- 


0,0+ 


0,0- 


1, 0+ 


i-,o 




14.16 


11.7 


0.43 


7.4 


2.0 


1.67 


_ 


coupling 


PV 


V 


V 


S 


PV 


S 


V 



Table 1: The masses M, the isospin 7, spin and parity J'^, and couplings g of non-strange mesons. The 
ratio of the tensor to vector coupling f /g — ior the cu meson and 5.1 for the p meson. 

meson fields lead to the following set of Schwinger- Dyson equations for nucleons 

G{p) = Go{p) + Go{pMp)G{p), (93) 
S(p) = -igl j -^,l,D>^'{q)G{p-q)T,{p-q,p-q) (94) 

- iglJ-^,Ds{q)G{p-q)r{p-q,p;q), (95) 

j^,^'^,TWDpv{q)G{p - q)rl\p - q,p; q), (96) 

where Go{p) and G{p) are the free and full nucleon propagators. The summation over the mesons 
with the same type of coupling is implicit. The nucleon self-energy S(p) contains the vector, scalar, 
and pseudo- vector meson propagators and the associated three-point meson nucleon vertices T{p). The 
meson propagators obey the following Schwinger-Dyson equations, which we write explicitly for the 
case of vector coupling, 

D'^'^iq) = D^''iq)+D^\q)U'^{q)D^^ (97) 

U^'-iq) = -igl j ^Jr[YG{p)Y{p,P + q;q)G{p + q)]. (98) 

The vertices T{p) obey their own Schwinger-Dyson equations which connect the three-point functions 
to four-point and higher order functions. The lowest order truncations of this hierarchy (i. e. replacing 
the vertices T{p) by their bare counterparts) leads to the rclativistic Hartree and Hartree-Fock theories. 
Another common approximation is to replace the meson propagators by their free-space counterparts; 
the resulting nucleon self-energy is written as 

" I W)~' " Ds{q)G{p - q)] 

- iglv J -0^,lVlVq''q''Dpv{q)G{p - q). (99) 
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Note that pions (which couple by the pseudo-vector couphng) contribute to the self-energy only via 
the Fock exchange term in last line of Eq. (j99|) . One recovers the conventional relativistic mean field 
models upon dropping the Fock exchange terms. (It should be noted that the meson self-interactions, 
which we neglected from the outset, play an important role in the relativistic mean- fields models. 
The self- interaction coupling provides a further tool for adjusting the models to the phenomenology). 
The phenomenological models that are based on the Hartree (or Hartree-Fock) description of nuclear 
matter (the theory is known also as quantum hadrodynamcis) have been used extensively to study the 
properties of nuclear matter; we will not discuss these models here (see the monographs (11 ^ .12; ^ .92 ^ .9,1]). 



2.4.2 Covariant T-matrix 

The theories which are based on the covariant treatment of the T-matrix and self-energy in nuclear 
matter are know as the Dirac-Bruckner-Hartree-Fock (DBHF) theories. These theories were developed 
during the last two decades mostly in the zero-temperature and quasiparticle limits jHEl UHl HUUl IIOH 
una [ID3 CDll CDSl IM IMl una una inni- This section gives a brief overview of the ideas underlying 
this theory. 

We start our discussion of the relativistic T-matrix theory by writing down the four- dimension 
Bethe-Salpeter equation (BSE) in the free space 

T(piP2;p[p'2) = V{pm-M) + ^ / ^^(piP2;p:,p'!)Gz5(p;)G'^(p':)t(p;,p'!;p;p'2)> (lOO) 

where P = pi +p2 is the center of mass momentum and p± = P/2±p, and Gd is the Dirac propagator 

iGoij)) = [y ■ p — m + iri]~^ (101) 

The reduction of the four- dimensional BSE to the three-dimensional form requires certain constraints on 
the zero-components of the four-momenta of the in- and outgoing particles. These constraints, the first 
one due to Gross [HI] and the second one due to Logunov-Tavkhelidze jHSI, Blankenbecler- Sugar 
and Thompson ^7\, require that 



p', = Jp'i + m^, pl = yG-pl (102) 



P? = P^ = ^, (103) 

where s = {pi +^2)^ is one of the Mandelstam invariants. The three-dimensional reduction of the BSE 
within the Thompson prescription is written as 

j^,V{p,p2-,pl,p'L)G2{p'',s)Tipl,p'L-,p[p'2-,s), (104) 

where for P = the two-particle propagator is 

G,(p..)^-|! ^;w-^;'-p) . (105) 

y/s - 2Ep + IT] 

where Ep is the on-shell particle energy and A^(p) are the projectors on the positive energy states 
(the negative energy states are commonly neglected, although a complete analysis of the covariant form 



30 



of tlie nucleon-nucleon amplitude requires information for both positive- and negative-energy Dirac 
spinors |,98j). 

After the Thompson reduction the interaction is instantaneous, i.e. the retardation effects intrinsic 
to the full BSE are removed. The reduced relativistic scattering two-body problem is thus described 
by BSE ()104|) which permits one to adjust the parameters of the interaction to the experimental phase 
shifts; the bound state spectrum is described by the homogeneous counterpart of Eq. ()104|) and can be 
used to constrain the interactions to reproduce the deuteron binding energy. 

Now we turn to the scattering problem in nuclear matter and write the formal solution of the 
Schwinger-Dyson equation for nucleons as 

zG(p) = [7-p-m- S(p)]~\ (106) 

The self-energy has a decomposition in terms of the Lorentz invariants 

S(p) = E5(p) + l.Kip) + ^S^^(p) + 75SP5(P) + l5l,KviP) 

= S5(po,p') + 7MP'^Sy(po,p') + ^p^p-1:t{po,p'). (107) 

The assumption that the theory is invariant under parity transformations requires that the terms 
involving 75 vanish; the last term in the second line vanishes by the anti-symmetry of the tensor a^^- 
Since we are interested in the equilibrium properties of matter, we shall not carry along the Schwinger- 
Keldysh structure and will specify the discussion to the retarded propagators. Upon separating the 
zero-component of the vector self-energy, 7^p^Sy(po,P^) = ~1o^o{p) +7'PSy(p), the propagator (|106|) 
can be written as a quasi-free (retarded) propagator 

iG{p) = - m*{p)]'^ , (108) 

where the effective momenta and masses are defined as 

m*{p) = m + J:s{p), pI = Po + Mp), P* = p[1 + ^v{p)] ; (109) 

the analogy to the Dirac propagator is formal because the new quantities are coupled via the self- 
energies and are complex in general. The form of the new propagator (jl08|) suggests defining effective 
spinors which are the on-shell positive energy solutions of the medium modified Dirac equation where 
the imaginary part is set to zero, i. e. 

where Xr is a state-vector in the spin-space and = ■>/ p*^ + m*^ is the energy eigenvalue. The effective 
spinors are normahzed according to Ur{p)Us{p) = Ssr- For further purposes it is useful to define the 
effective quantities according to E{p) = {m* / E*)u{p)E{p)u{p). Acting on equation ()108p by the unity 
operator A"*"* + A~* = 1, where the positive and negative energy projectors are defined as A"*"* = Ur®Ur 
and upon neglecting the negative energy part one finds 

G{p) = [pl-E; + ^C{p)y\ (111) 

where the damping is defined as 



C{p) = Im 



^o{p) - T^^sip) - l^pSy(p) 



(112) 
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The spectral function can be constructed in full analogy to the non-relativistic case 



The second relation, which corresponds to the quasiparticle limit of the spectral function, defines the 
single particle energies 

£p= /FFT^-So(p). (114) 

As in the non-relativistic case, the form of the spectral function is Lorentzian and the spectral sum rule 
(j^^ is fulfilled for the general and quasiparticle forms of the spectral function in Eq. pi3|) . 

The Bethe-Salpeter equation in the background medium and in the reference frame of the center- 
of-mass of two particles (suppressing spin indices) is written as 

TM;s',n - VM) + / ^Vi,, ^ 1'-^^' ^ „ n.,v'; s', P) (115) 

where Q2{(1, s*,P) = 1 — fF{£p/2+q) — fFi^p/2-q) is the Pauli blocking and ± is a short hand for P/2±p. 
The dependence of the Pauli blocking on s* and P is due to the function evaluated in the two-particle 
center-of-mass frame, where the Fermi-sphere is deformed because of Lorentz transformation from the 
lab to the center of mass frame. A closed set of equations is obtained upon introducing the retarded 
self-energy in terms of the T-matrix (jll5|) 

s(p) = I ^/rD[Tr(r + m'^)/(2)]T° - (r + mnf^2)T:j, im 

where the subscript ex stands for exchange, T" are coefficients of the expansion of the full T-matrix in 
Lorentz invariants 

^ = E /rD/(2)^^ f(^) e {1, l^ll 75 4^]. (117) 

a 

The chemical potential appearing in the Fermi-functions is adjusted to reproduce the density of the 
system. The solutions of the self-consistent, finite-temperature relativistic T-matrix theory allow one 
to compute the energy density as 



E 



(2^ 



{u*{p)\-fp + m+^j:{p)\u*{p)) 



fpisp), (118) 



and the thermodynamic quantities introduced in Eqs. (j42j) and ()44j) . The binding energy at zero- 
temperature [fp{ep) = 9{pf — p)] is obtained from Eq. (|118j) 

Eb{pf) = P~^E{pf) - rriN. (119) 

While we have kept only the positive energy states in our discussion, an unambiguous treatment of 
the nucleon self-energy in matter requires keeping the negative energy states as well llU5t llU6j . An 
example of such an ambiguity is the pion exchange part of the Lagrangian which can be described by a 
pseudo-scalar or a pseudo- vector coupling. Both couplings produce the same free space matrix elements 
for the on-shell nucleons when the coupling constants fpy and gpv are related as fpv/gps = 171^, /2m, 
where is the pion mass. If only the positive energy states are kept, a recipe to overcome this problem 
is to divide the T-matrix into the Born term plus a correlation term |107| I108j . The ansatz ()117p is 
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Figure 13: Left panel: The non- vanishing components of the self-energy S5, Sq and 
as a function of the Fermi-momentum within the Dirac-Brueckner (heavy hues) and the 
relativistic T-matrix (hght hues) approximations. Right panel: The binding energy per 
particle of isospin symmetrical nuclear matter within the Dirac-Brueckner and T-matrix 
theories |lU2j . The dashed line shows binding energy in a non- relativistic Brueckner theory 
where the three-body forces are neglected. The empirical saturation point is shown with 
error bars. 



applied only to the correlation term, since the structure of the Born term is dictated by the interaction 
V, which is fixed. 

Another often used approximation is the neglect of the momentum dependence of the self-energies, 
which are approximated by their value at the Fermi- momentum pp. The form of Pauli-operator in Eq. 
flllSj) which keeps only the particle-particle propagation in the intermediate state is the counterpart 
of the non-relativistic Bruckner theory (as it relies on the ideas of the hole-line expansion) jlOm I1UX| . 
Including the hole-hole propagation leads to the relativistic counterpart of the original T-matrix the- 
ory [Tn2j . 

Fig. El left panel, shows various components of the nucleon self-energy from a simplest type cal- 
culation which ignores the momentum dependence of the self-energies, the negative energy sea and 
works at zero temperature |lU2j . The contribution of Eq component is negligible in the case where the 
negative energy contributions are neglected and the nucleon-nucleon amplitude is expanded according 
to the ansatz ()117j) . The components S5 and Eq are large on the nuclear scale, but of the same order 
of magnitude, so that their contributions mutually cancel. The binding energy of isospin symmetrical 
matter with the Bruckner and T-matrix approaches is shown in Fig E[ right panel. The additional 
density dependence of the Dirac spinors is an important ingredient of the relativistic T-matrix theories 
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which leads to a new saturation mechanism. Compared to the non-relativistic theory the saturation 
density is correctly reproduced by the relativistic theories. The role of the three-body forces within 
the relativistic theories, in particular the role played by the A isobar, has not been addressed in the 
literature. 

2.5 Isospin asymmetric matter 

The proton fraction Yp in neutron star interiors is constrained by the condition of equilibrium with 
respect to the weak processes. The disparity between the neutron and proton numbers (breaking the 
SU{2) symmetry in matter) motivates the study of the nuclear matter under isospin asymmetry, which 
conventionally is described by the asymmetry (or neutron excess) parameter a = (p„ — Pp)/{pn + Pp), 
where pn and pp are the number densities of neutrons and protons, or alternatively by the proton fraction 
Yp = (1— a)/2. The isospin asymmetry is accommodated in the T-matrix and related theories by working 
with two-point functions (self-energies, T-matrices etc.) which are 2x2 matrices in the isospin space. 
The fundamental quantity characterizing the asymmetric nuclear matter is the symmetry energy, Es 
(i. e. the energy cost of converting a proton into a neutron). For small values of a the symmetry energy 
can be expanded in series 

Es = S2a'^ + S^a^ + 0{a^), (120) 

where 5*2 is the coefficient in the symmetry energy term of the Bethe-Weizsacker formula. At zero 
temperature the contribution from the kinetic energy to Eg can be evaluated explicitly in the theories 
where the matter effects are included in interaction energy alone 

Ej^ip, a) - Ej^ip, 0) = EjskS.^ [(1 _ «)5/3 + (1 + ^)5/3] _ i| ^ 5^^^^^ ^^^2 ^ ^^^4). (^21) 

the coefficient of identifies the contribution of the kinetic energy to 5*2. The interaction energy con- 
tribution is clearly model dependent, in particular it depends on the many-body theory and the nuclear 
interaction adopted. The importance of studying the symmetry energy arises from the importance 
of neutron /3-decay reactions in high density matter, whose density threshold depends on the proton 
concentration (see SubsecE^- Some constraints on the symmetry energy can be obtained around the 
saturation density po- An expansion of 5*2 with respect to small deviations from po gives 

S,{n) = SM + ^(P - Po) + 0[{p - pon (122) 
Po 

The first derivative of 5*2 determines the change in the pressure at the saturation point due to the 
asymmetry of the system. Upon using the expansion p2()|l one finds ['T^^ 

= . P.. (123) 

The symmetry energy as a function of density is shown in Fig. for several models which are based on 
the relativistic Dirac-Bruckner approach 111211113] . variational approach |114j and BBG approach |115j : 
the latter two approaches include the three-body forces. 

The values of the symmetry energy and its derivative at the saturation density vary in a narrow 
range: S'2(po) ~ 29 ± 2 MeV and Pq = 3 ± 1 MeV fm^. The predictions of various models of the high 
density behavour of the symmetry energy differ substantially. The relativistic, variational and BBG 
theories (the latter without three-body force) vary within 10% of an "average" value. The BBG theories 
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Figure 14: Dependence of symmetry energy of nuclear matter on density. CT[1] and CT[2] 
refer to the results of refs. |112j and fL13\ obtained within the covariant T-matrix theory. 
V refers to the variational result of ref. |114j which includes a three-body force; BBG refers 
to the results based on the BBG theory with two-body forces |115j : BBG(l) includes in 
addition a microscopically derived three-body force IllTl I118j : BBG (2) is the same as 
BBG(l) but with a phenomenological three-body force |1191 112U] . 



supplemented by either microscopic |116| I117t I118j or phenomenological |119t ll2Uj three-body forces 
predict symmetry energy that is by a factor of two larger than predictions of other models. However, 
the discrepancies in the magnitude of the symmetry energy at asymptotically large densities are not 
essential, since other degrees of freedom such as hyperons, mesonic condensates, or other states of 
matter are likely to occupy the stable ground state. 



2.6 Hyperons 

At densities around the saturation density the only baryonic degrees of freedom are protons and neu- 
trons, which form an iso-duplet whose approximate free-space SU{2) symmetry is largely broken in 
matter. At larger densities the number of stable baryons increases. These include the isospin 3/2 nu- 
cleon resonances A^, A° and the strangeness carrying baryons (hyperons). The hyperonic states can be 
classified according to the irreducible representation of the SU (3) group. The two diagonal generators 
of the group are linear combinations of the isospin T3 and hypercharge Y, which is equal to the sum 
of the baryonic number and strangeness, Y = B + S. The charges of baryons are determined by the 
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Figure 15: The SU{3) baryon octet. The axis are the isospin T3 and the hypercharge Y. 
Baryon masses are given in brackets. 



Gell-Mann-Nishijima formula Q = Y/2 + T^. Fig. 1151 shows the octet of the baryons whose strangeness 
carrying members can appear in neutron star matter. If we neglect the interactions between the hy- 
perons and nucleons, the threshold for hyperons to become stable is determined by comparison of the 
hyperon mass to the largest available energy scale - the neutron Fermi energy |121j . The hyperons 
can appear in matter through the weak hyperonic (inverse) beta-decay reactions e~ + n + z/g and 

hadronic weak decay n + 7r~ — > S~. The energy balance in the first reaction implies 2/i„ ~ M^- = 1197 
MeV, where /i„ is the chemical potentials of neutrons (we used the fact that the chemical potential 
of neutrons and electrons are almost equal in matter under (5 equilibrium, see subsection 12. 7^ . The 
r. h. side of the second reaction is 0{fin), therefore it is negligible compared to the first reaction. 
Similar arguments apply to other hyperons which are stabilized either through the hyperonic [3 decays 
or hadronic weak decays. For example for the lightest hyperon A" one finds 

p + e- - AO + u p + n-^ A^, 0(^,) ] ^ ^ 1116 MeV (124) 

The reactions in the first line being 0{fip), where /ip <C /i„ is the proton chemical potential, can be 
neglected and A^ appear primarily through the weak hadronic process in the second line. Since the 
r. h. side of this reaction is 0(/i„) and the mass difference Ms- — Ma < //„ at relevant densities, S~ 
hyperons appear first. 

Interacting hypernuclear matter was initially studied within variational approaches by Pandhari- 
pande |131j and Bethe and Johnson |132j . With the advent of the powerful phenomenology of relativistic 
mean-field models these were extended to the hyperonic sector |133| 1134(11351 1136] . The extension of the 
T-matrix and related theories to include hyperons requires the interactions between the hyperons and 
nucleons (YN) and hyperons and hyperons (YY) \1'2'2\ 11231 11241 H"^ - The experimental information 
about the interaction involving hyperons is rather scarce. The YN potentials are fitted to the AA^ and 
SA^ scattering data. The information on YY interactions is limited to the ground state of double-A 
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Table 2: The interaction channels within the isospin basis, for combinations of the total strangeness S 
and total isospin /. 

7 = / = i 1 = 1 1=1 1 = 2 



S = NN NN 

s = -i {an,j:n) j:n 

S=-2 (AA,SiV,SS) (SiV,SA,SS) ES 

S=-3 (SA,SS) SS 



hypernuclei |126j . Additional constraints come from the SU{3) symmetry arguments. The main dif- 
ference between the YN and the ordinary NN interactions is that the direct YN interaction does not 
contain the one-pion-exchange (hereafter OPE) part of the A^A^ interaction, therefore the short range 
part of the nuclear force is not hidden under the dominant OPE interaction. The A hyperon couples to 
the neutral pion due to A — S° mixing. 

The extension of the T-matrix theory to include the hyperonic degrees of freedom requires a treat- 
ment of the coupled-channel problem. The possible interaction channels in the isospin basis are given 
in Table |21 In the strangeness S = 0, —4 sectors there is a single channel. In the S = —1 sector the 
channels AA^, T^N are coupled and the T-matrix equation reads 

TaN;AN TaN;T:N \ _ f Van^AN VaN;T:N 

VaN;AN VaN;T,N \ ( GanTaN:AN GanTan:T,n 



where the intermediate state propagator, which generalizes the single species result fl33|) to a multi- 
component system, is 

i2^' jj^^MpM-) II - MP.) - /B.(P-)| J;^, (126) 

where B stands for any baryon of the S'f/(3) octet. The spectral functions in Eq. ()126|1 are related to 
the self-energies 

AN T,N ' ' ^ ^ ^ ' -'A^^'■A^^' -'■A^r■v^r \ / i\ I '^\at >-J 



= E T^^''''' T^^'"''' ]{P + P)[ ]{P') 

+ 5^2,(u; + u;')f + p')Im f Vp')- (127) 

Equations (jl25j) and (jl27p are the generalization of the equations of (jH^ and (PT|) to the case of coupled 
AA^, SA^ channels. The example above is sufficiently general to illustrate the treatment of other coupled 
channels shown in Table El The Brueckner-Bethe-Goldstone theory for hyperonic matter is recovered 
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from the above equations by (i) taking the quasiparticle and zero temperature hmits and (ii) dropping 
the hole-hole propagation from the intermediate state propagators. 

The BBG calculations of hypernuclear matter were carried out over several decades in parallel to 
the development of the theory for ordinary nuclear matter with a special attention to the problem of 
binding of A particle in nuclear matter (see e.g. |127t 11281 H^i^l H^Oj ^. Recent work on Bruckner theory 
for hypernuclear matter shifted the interest towards understanding the /3-equilibrated matter in neutron 
stars dHZl CHHl Cnni En] (see SubsecEl below). 



2.7 Charge neutrality and weak equilibrium 

Neutron stars evolve towards equilibrium with respect to the weak interactions on long time scales. 
However, in many cases, the time-scales of interest are much shorter than these equilibration time- 
scale and we can assume, to a good approximation, that the NS interiors are in approximate weak 
equilibrium between hadrons and leptons. If the temperatures in the interiors of NS are below several 
MeV, which corresponds to time scales of the order of months or less after star's formation, the neutrinos 
propagate through the star without interactions. Their chemical potentials can be set to zero. In a 
matter composed of neutrons (n), protons (p) and electrons (e) the weak equilibrium is established by 
the /3-decay and electron capture reactions 

n— >p + e~ + z/e, p + e~ ^ n + Ue- (128) 

The equilibrium condition requires that the chemical potentials obey the equality = fJ'n — f^p- Given 
the energy density of asymmetric nuclear matter e the chemical potentials can be expressed as 

de de 1^ a de ^ 2 de . . 

oun/p on n oa n oa 

where 6fj, = fj,n — fip = fie- The charge neutrality requires that the number densities Up = Ue] treating 
electrons as a non interacting and massless rie = /ig/Svr^, which in turn implies that 

9eV 8 os/l 



where in the last step the expansion p20|) has been used. Eq. ()130|) allows one to compute the proton 
fraction at a given density if the symmetry coefficient 5*2 is known. When the electron chemical potential 
exceeds the muon rest mass = 105.7 MeV the electrons decay into muons via the reaction 

^ fi' + + Ue- (131) 

Equilibrium with respect to this process and its inverse requires /Xg = /^^t, where fi^ is the chemical 
potential of muons. The proton fraction in this case is given by the parametric equation 

y, = J^[(5^^-my3/2 + 5^3]. (132) 

The conditions for /5-equilibrium and charge neutrality are readily generalized to the case of arbitrary 
numbers of baryons and leptons; the key observation is that there are two conserved charges available 
- the total baryonic charge, which is related to the conservation of the baryonic density and the total 
electrical charge, which is related to the charge neutrality of matter. The thermodynamic potential 
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Figure 16: Matter composition in /9-equilibrated hypernuclear matter. Left panel: Noninter- 
acting matter |121j : Central panel: Interacting mixture according to the BBG theory |14Uj : 
Right panel: Interacting mixture within a relativistic mean- field theory |133j . 



of matter is a functional of the baryon n{Yi) and lepton n{Yj) densities, where i and j enumerate the 
baryon and lepton species respectively: 



n{Y„Y,) = E - fXB 



riB 



Y^BiY.n 



(133) 



where Bi and Qi are the baryon and electric charges and fiB and /i^ are the associated Lagrange 
multipliers, is the baryon density. The equilibrium conditions require 



Si', • dY, 

The baryon and lepton chemical potentials are obtained as 



dE o ^ 



dE 
dY 



(134) 



(135) 



If the neutron and electron chemical potentials are chosen as independent parameters, Eqs. (jl35|) 
identify /is = —fin and fie = —fJ^L and the chemical potentials of arbitrary baryonic species are written 
as lii = Bijin — Qijie- For example for a mixture consisting of n, p, S='= and A° baryons the equilibrium 
conditions give 

Yp + ys+ - {Ye + Y^ + Yj,-)= 0. (137) 



39 



Fig. El shows the abundances of various species in a baryon-lepton mixture within the non-interacting 
hyperonic matter theory |121j . BBG theory jl4Uj and the relativistic mean-field theory |133j . The 
interacting models agree quahtatively at low densities, however at larger densities the relativistic theory 
allows for the S^'"*" hyperons and the cascade S~. 



2.8 Meson condensation 

Bose-Einstein condensation (BEG) of mesons has an important impact on the properties of hadronic 
matter in NS. First, it softens the equation of state and, second, it leads to an enhanced neutrino 
emission. Contrary to the ordinary Bose-Einstein condensation, where above the critical temperature 
one deals with a normal Bose gas, the pion condensation is associated with an unstable mode in nuclear 
matter which has the quantum numbers of pions (ref. [141] and references therein). The threshold of 
pion condensation is derived by considering the pion retarded propagator in nuclear medium 

D^ioo, q) = [u^ -q^-ml- n(cu, q)] , (138) 

where the pion self-energy, which is represented by the polarization tensor of nuclear matter n(u;,g), 
sums the particle-hole {ph) and A-resonance - hole (A/i) states to all orders. [The equations determining 
the pion propagator and the self-energy are identical to Eqs. and (jHEI) with the vector coupling 
replaced by the pseudo-vector coupling.] The resummation of these channels in nuclear matter is a 
complex problem in general |142t I143j . but can be performed analytically if one approximates the 
driving term in the ph and Ah series by Landau parameters 



ph 



f + g(Ti-(T2+{f' + g'(Ti- (72) {Ti ■ T2), (139) 



where cr and r are the vectors of Pauli matrices in spin and isospin spaces. In analogy to p39p one 
may define the interaction in the Ah channel as Fa/i = A'ataI^'i ' S2)(ti • T2) and Faa = 5'aa('^i ' 
S2){Ti ■ T2), where S and T are the spin and isospin operators for the A resonance. The net 
polarization tensor is the sum of nucleon and A contributions, 11 = (u;^ — q^)xN + q'^Jxa, where 
/ = [{rriN + ttia)^ — (i^^ — q"^)] /4mA, with m^v and being the nucleon and A-resonance mass. The 
"susceptibilities" which include the effects of short-range correlations are given in terms of one-loop 
Lindhard functions Xn and xa as |144j 

Xn = [1 + {9'aa - 9'na)xa]V'^, Xn = [1 + {9nn - 9'na)xa]V'^, (140) 

where V = 1 + 9'^^XN + 9'aaXa + {9'nn9'aa + 9na)XnXa- Due to the isospin symmetry of nuclear 
matter the polarization tensors of neutral and charged pions are equivalent. Their spectral function can 
be written as 

-2Imn^(cj, q) 
~ - g2 _ ^2 _ Ren^(c^, g)]2 + [ImU^itu, g)]^ ' 

The pion condensation in symmetric nuclear matter is characterized by the condition B{0, q^) 00 
at lmn^(0, qc) ^ 0. In the presence of vr condensate the uniform nuclear matter acquires a periodic 
nucleonic spin wave structure with wavenumber q^. These qualitative features remain intact for neutral 
pion condensation in neutron matter. Charged pion condensation in neutron matter is characterized 
by additional instabilities. Apart from the collective modes mentioned above, a mode appears which 
carries the quantum numbers of vr"*" above some critical density p+, which depends on the details of the 
repulsive interaction in the S = 1 and T = 1 channel. In terms of nucleon excitations it represents a 
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bound state of a proton and a neutron hole. At higher densities the sum of the poles of the vr"^ and 
TT" propagators [c. f. Eq. (jl38|) ] vanishes, which signals the instability of matter towards formation of 
7r^7r+ meson pairs. 

Whether the pion condensation occurs in compact stars is an open issue. The answer depends 
crucially on the values of the g' parameters; the currently accepted range of g'j^j^ G 0.6 — 0.8 and the 
universality ansatz, which sets all the g' parameters equal to S'tvat; precludes pion condensation in finite 
nuclei and in compact stars within the density range where the nucleonic and mesonic components 
retain their identity. Recent analysis of Gamow- Teller resonance in the ^°Zr(p,n)^°Nb reaction [TIH] 
suggests the following values of parameters: g'j^j^ = 0.59 and g'j^^ = 0.18 + 0.05^^^^, with g'^^ being 
undetermined. With this new information the critical density of pion condensation turns out to be 
lower than that with the universality ansatz [146J. An independent evidence for low-density neutral 
pion condensation was obtained in variational calculations of ref. |114j . 

The mechanism by which kaons may from a Bose condensate in neutron star matter was developed 
by Kaplan and Nelson [HI]. Since the anti-kaon interactions in nuclear matter are attractive, their 
effective mass could be substantially lower than their mass in the vacuum. Thus, instead of neutralizing 
the positive charge by the negative charge of energetic electrons, this can be done by stabilizing K~ 
in matter. Most of the effects of K~ condensation on the properties of compact stars are similar to 
those discussed for pion condensation. These include a substantial softening of the equation of state 
which reduces the maximum mass of a compact star to 1.5Mq. This reduction may have important 
implications for the low-mass black hole population in our Galaxy |147j . 



2.9 Stellar models 

The equilibrium configurations of compact stars are described by the Einstein equations of General 
Theory of Relativity (GTR) 

1 Sirk 

Rik - -Rgik = ^Tik 2, A; = 0,1, 2, 3 (142) 

where Rik is the Ricci tensor, R - the scalar curvature, gik - the metric tensor, Tik - the energy-momentum 
tensor and k - the Newton's constant; the cosmological constant is omitted from the Einstein's equa- 
tions. The gravitational field created by a spherically symmetric mass distribution is itself spherically 
symmetric. For such fields the components of the metric tensor are functions of the radial coordinate 
r and time t. The metric can be written 

ds^ = e'c'dt^ - e^dr^ - r''{de^ + sin ''Odcf)''). (143) 

The energy and momentum tensor can be expressed through the mass density p and pressure P as 

T,fc = {p + c-^P)uiUk - Pgik, (144) 

where Ui is the matter four-velocity. For a static mass distribution Tq = pc^ and = = T| = —P. 
In the static (time- independent) limit and for spherically symmetric gravity Eq. (jl42j) reduces to 

^ {---^]+-^ = -TP (145) 
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where primes denote radial derivatives, 
dynamic equilibrium 



The last equation can be replaced by the equation for hydro- 



0, 



(148) 



which is nothing else than the explicit form of the covariant hydrodynamic equations = 0. Eqs. 
()145j) . p4(i|l . and p48j) should be supplemented with the equation of state -P(p) and boundary conditions 
P(0), z/(0) and A(0). The inner solution within the configuration should be matched with the external 
Schwarzschild solution 

= e-^ = 1 - 149 

where M is the mass of the configuration. The quantity Rs = 2kM/c^ is the Schwarzschild radius; 
the importance of general relativity is determined by the ratio Rs/R where i? is a characteristic length 
(e.g. the star radius). 

It is evident that at the center of the configuration A(0) = 0, A'(0) = z^'(O) = to avoid singularities 
as r 0. The problem of finding the internal solutions of stellar configurations is simplified upon 
introducing a new variable via the relation = 1 — 2km{r)/c^r. One finds |148j 



dm{r) 

dr 
dP 

dr 



k{P + cV) 
[c^r — 2km{r)]r 



m{r) H — :rPr^ 



(150) 
(151) 



Thus we have two equations for three variable p, p, and m. The system of equations is closed by specify- 
ing the equation of state P{p). We implicitly assumed that the pressure is independent of temperature, 
as is the case for a system at zero temperature. However in those cases where the temperature is impor- 
tant, the complete set of equations includes apart from the equation of state also two new differential 
equations which describe the energy balance and thermal transport equations (see Chapt. 1^^ . 

To gain insight in the equations of hydrostatic equilibrium ()15()j) and ()151|) consider a self-gravitating 
star with constant density po- Integrating the first equation gives m(r) = 47rpo'"^/3. The integration of 
the second equation gives 



Rq 



1 - 



3P + cVo Pq + cVo 
P + cVo 3Po -t- c2po 



1/2 



Rn 



3c' 



Srrkpc 



1/2 



(152) 



Substituting P 
as 



one finds the radius R of the configuration. The metric within the sphere is written 

2kM\ 



^ r 



1 - 



cm 



■ 



r < R. 



(153) 



^c^po + P(r 

These solutions match the Schwarzschild solution ()149j) at the surface of the configuration. For the 
constant density po = 2.85 x 10^^ g cm~'^ and maximal pressure P = 5 x 10'^^ erg cm~^ one finds the 
maximal radius and the mass of the configuration: R = 6.4 km, M = 0.16 Mq. In the case of an 
incompressible fluid P{r)/c'pQ = 1/3 and one finds R = 17.74 km, M = 3.3 Mq. The numbers above 
give approximate lower and upper limits on the mass and radius of a compact object. For a more 
realistic estimate we now turn to the many-body equations of state of hadronic matter. 

For illustrative purposes we have chosen several equations of states which are either based on the 
theories discussed in previous sections or are alternatives listed in Subsec. 12.101 their properties are 



42 



5e+14 le+15 1.5e+15 2e+15 2.5e+15 3e+15 

% [gem ] 

Figure 17: Pressure vs central energy density for representative equations of state (EoS). 
Model A [solid line) - relativistic mean field model, features hyperons and pion conden- 
sate ■ Model B [dashed line) - non-relativistic variational EoS based on AV14 potential 
supplemented by UVII three-body force |149j . Model C [short dashed line) - relativistic 
model based on covariant T-matrix |112j . Model D [dashed double- dotted) - non-relativistic 
EoS based on the BBG theory, includes a three-body force. Model E [dashed dotted) - same 
as model D, but includes interacting hyperons |14Uj . 



summarized in Table |21 Dependence of pressure on central density for models labeled A-E are shown 
in Fig. El Model A is a relativistic mean field model which allows for hyperons and onset of pion 
condensation |133j . Model B is a non-relativistic variational model based on the Argonne AV14 two- 
body potential which is supplemented by the Urbana UVII three-body force |149j . Model C is a 
covariant model based on the Dirac-Bruckner (covariant T-matrix) theory |112j . Finally, models D and 
E are based on the BBG theory which employs two-body and three body forces |14Uj . Model E includes 
hyperons along with the nucleonic component, whose interactions are taken into account within the 
BBG theory 

Equations of state are commonly characterized by their stiffness, which is expressed through the 
adiabatic index 7 = logP/logp (the larger the adiabatic the stiffer the equation of state). The non- 
relativistic models B and D which include only nucleonic component interacting with two- and three- 
body forces are the stiffest equations of state at large densities. The relativistic mean field model A 
is stiff at low densities where the hyperons are still not present but becomes softer at higher densities 
as the result of hyperonization of matter and pion condensation. As a general trend new degrees of 
freedom soften the equation of state since the new constituents of matter share the stress due to the 
pressure with other constituents. Such a softening is apparent from a comparison of non-relativistic 
models D and E which differ by the presence of hyperons in the latter model. Note that the models 
are shown in their published density range. They were supplemented by the equations of state for 



43 



EoS method composition forces 



A 


RMF 


npeHn 


contact 2E 


B 


variational 


npe 


realistic 2E + 3E 


C 


DBHF 


npe 


realistic 2E 


D 


BEG 


npe 


realistic 2E + 3E 


E 


BEG 


npeH 


realistic 2E + 3E 



Table 3: Summary of EoS A-E introduced in the text; H refers to hyperons, vr - pion condensate, 2B 
and 35 - two and three-body forces, RMF - relativistic mean field. The remaining abbreviation are 
introduced in the text. 



EoS 




M 


R 




Pk 


T/W 


cJ/kMl 


/ 




[10^^ g/cm^] 


Mq 


[km] 


[10^ s-1] 


[10-3 s] 






[1045 g cm^] 


A 


4.00 


1.98 


13.1 


1.05 


0.60 


0.088 


2.32 


1.94 


E 


2.00 


2.50 


12.9 


1.20 


0.52 


0.145 


4.53 


3.29 


C 


2.00 


1.77 


12.9 


1.03 


0.61 


0.118 


2.08 


1.77 


D 


1.60 


2.49 


13.5 


1.14 


0.55 


0.181 


4.67 


3.62 


E 


1.60 


1.47 


12.8 


0.97 


0.64 


0.200 


1.55 


1.41 



Table 4: Parameters of compact stars rotating at their Keplerian frequency constructed from EoS A-E 
specified in the text. The models listed are the largest mass (stable) configuration for the EoS A and E 
and the largest central density configuration for EoS C-D. The entries are the central energy density, the 
gravitational mass, the equatorial radius, the Keplerian rotation frequency and period Pk = 27r/f2j^, 
the ratio of the kinetic (rotational) to the potential (gravitational) energy, the angular momentum J, 
and the moment of inertia /. 

low-density matter from refs. |150j and |151j to construct sequences of stellar models. 

The masses and radii for sequences of configurations with different central energy densities Sc are 
shown in Fig. for equations of state (hereafter EoS) A to E. The configurations were computed with 
the RNS code written by Stergioulas |152j . As a common trend one finds that the stiffer is the EoS 
at the central density of a configuration, the larger is its mass. Each configuration features a maximal 
mass, but because of the limited density range across which our models are defined this mass is apparent 
only for models A, E and D . Mass measurements of compact stars obtained from the timing of binary 
pulsars are broadly consistent with a canonical mass of 1.4 Mq. Pulsars that have undergone long 
periods of accretion have larger masses. The most massive measured pulsar to date PSR J0751+1807 
is a millisecond pulsar in a 6hr binary system with a helium white dwarf secondary. The mass of the 
pulsar is measured through the orbital decay which is interpreted as due to emission of gravity waves. 
Combined with measurements of the Shapiro delay, this implies a pulsar mass 2.1 ±0.2 Mq |153j . While 
models A-D feature stars with canonical masses of 1.4 Mq (an exception is the model E), the latter 
observation places a severe limit on the EoS of nuclear matter and excludes all the models except E and 
D. Recent observations of binary pulsars that have undergone extended periods of accretion without 
losing their stability indicate that the EoS of compact stars ought to be rather stiff. All except the 
lightest members of the sequences shown in Fig. are characterized by small radii which are confined 
within the range of 8-11 km, their value being almost independent of the central energy density (see 
Fig. IT5|) . This leads to degeneracy in the mass-radius relationship in Fig. for models E, D and 
E there exist configurations with different masses but same radii. Apart from the stationary global 
observable like the mass and radius, stellar configurations are characterized by the dynamical quantities 
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Figure 18: Left panel: Dependence of the gravitational mass (in units of solar mass Mq) 
on the central energy density for configurations constructed from EoS A-E specified in the 
text. Right panel: Dependence of the star's radius on its central density for EoS A-E. 



such as the moment of inertia J, rotational angular momentum J and limiting rotation frequency ^Ik, 
know as the Kepler frequency. The rotation of the star modifies the metric and thus the rotation 
frequency in the local inertial frame, which we denote by ujl- Slowly rotating stars admit perturbative 
approach, where the small parameter is the ratio of the rotational kinetic energy to the gravitational 
binding energy |154t I155t I156| I157| I158t I159j . The angular velocity of a slowly rotating star obeys the 
following equation 

(r^jcj^)' + Ar^fuL = 0, j = 3'^"+^^^^ (154) 
The external solution of Eq. ()154j) is given by 

u^,{r) = - ^ (155) 

where Vl is the rotation frequency. The internal solution can be obtained through integration of Eq. 
()154p with the boundary conditions tu^ = 0, when r = and j{r) = 1 when r = R. The moment of 
inertia for slowly rotating objects can be written as 

^ ^\rr'e^'-^y'(p+^]^, (156) 



where the volume integration assumes approximate spherical symmetry. The Keplerian frequency Qk 
is defined as the frequency at which the centrifugal force on a test particle at the equator equals the 
gravitational binding force. For rotational frequencies Q > Qk mass shedding from the equator leads 
to instability of the configuration. Tabled lists the parameters of the general relativistic configurations 
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Figure 19: Dependence of configuration masses on their radii for EoS A-E specified in the 
text. 



rotating at the Kepler frequency. The values of the central density correspond either to the maximum 
mass stable (non-rotating) configuration (models A, B) or the largest central density available for a given 
EoS (models C-E). N. B. The central densities of these configurations exceed the nuclear saturation 
density by factors of the order of 10. It is likely that the baryonic and mesonic degrees of freedom are 
not appropriate at such large densities; nevertheless the EoS of quark matter can not differ dramatically 
from those listed above, and the conclusions drawn from the limiting frequencies should hold, at least 
qualitatively. The fastest millisecond pulsar measured to date is IGR J0029 1+5934 with rotation period 
P = 1.7 ms. According to Tableland the studies of limiting frequencies of rotating superdense stars 
based on alternative EoS (ref. ^2); Chapter 16), stars that are gravitationally bound can not rotate 
faster than half a millisecond. The situation is different if a configuration is self-bound due to strong 
interactions, as is the case for the strange stars. Thus, if an object is observed in the future with 
rotation periods smaller than half a millisecond it must be an exotic object, e.g. a compact star made 
up of strange matter (ref. [12], Chapter 18). 

2.10 A guide to alternative methods 

Here we complement the discussion of the nuclear many-body problem above by a brief summary of 
some of the alternative methods that are used in the studies of the nuclear many-body problem. 

Qunatum Monte Carlo (QMC) mei/iO(is have been used extensively to study the properties of infinite 
nuclear matter, neutron matter and light nuclei with various versions of the Argonne two-body and 
the Urbana three-body interactions (ref. |16(Jj and references therein). Recent developments include 
calculations of the ground state energy of nuclear and neutron matter using the Argonne AV18 and 
Urbana UIX potentials |114j . Green's functions Monte-Carlo (GFMC) calculations of light nuclei |161j . 
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and QMC calculations of neutron matter |162l 11631 1164j . These calculations are variational in nature 
(i. e. provide an upper bound on the energy of the system) and are based on a non-relativistic description 
of nucleonic matter in terms of the Schrodinger equation. 

Chiral Perturbation theory seeks to establish a connection between non-perturbative QCD and low- 
energy phenomenology of nuclear systems. This top to bottom (with regard to the energy scales) 
approach is anchored in the symmetries of the QCD Lagrangian (notably the approximate chiral sym- 
metry of strong interactions) and in the QCD sum rules. In the nuclear and neutron matter problems the 
chiral perturbation theory offers a method of treating the long-range pion-nucleon dynamics explicitly, 
the short-range correlations being incorporated in contact terms |165| ll6bT I167t ITBHj . As mentioned in 
Subsec. 12.81 chiral Lagranians are useful in deriving the properties of the kaon condensate in dense mat- 
ter. Furthermore, chiral Lagrangians have been used to derive free-space nucleon-nucleon interactions 
that are input in many-body and few-body calculations (ref. |169j and references therein). 

Relativistic density functional theories are based on the ideas of the mean-field theory of nucleons 
and mesons. These models incorporate a limited number of phenomenological constants, that are fitted 
to the properties of bulk nuclear matter and finite nuclei, and provide a powerful tool to study many 
aspects of nuclear phenomenology at an elementary (Hartree or Hartree-Fock) level |17Uj . Some recent 
models incorporate the chiral symmetry in the Lagrangian of the theory [1711 1172[ I173j . 

Lattice field theory methods have gained attention in recent years. The nuclear and neutron matter 
problems were studied in close analogy to the numerical simulations of the Hubbard model in condensed 
matter systems |174t I175| 1176] . Lattice methods were also applied to dilute, spin 1/2 non-relativistic 
fermions |177j . A lattice realization of the scalar 0^ field theory was applied to study a clusters in 
nuclear matter |178j . 

During the past decade the Effective field theory methods were developed for nuclear systems and 
applied to the few-body nuclear physics; see refs. [182, 183, IS^j and references therein. Some aspects 
of the many-body problem are discusses within this method in refs. |185| I186j . 

Coupled cluster method is one of the most powerful and universally applicable techniques in quantum 
many-body theory, that is well suited for treating non-perturbative problems. Its application to the 
nuclear many-body problem and finite nuclei has a long history (see ref. [llYQj and references therein). 
Recent developments include the quantum chemistry inspired coupled cluster calculations of ground 
and excited states of nuclei |18()[ 1181] . 

3 Neutrino interactions in dense matter 

The temperature of a NS born in a supernova explosion is of the order of several tens of MeV. During 
the first seconds neutrinos are trapped inside the star (i.e. the neutrino mean free path <^ i?, where 
R is the star radius). The energy that is lost in neutrinos is radiated from the surface of the neutrino- 
sphere - the surface where the optical depth for neutrinos drops to zero. Once the star cools down 
to temperatures of the order of several MeV, the matter becomes transparent to neutrinos {\y ^ R). 
The subsequent thermal evolution of NS is controlled by neutrino radiation from its interiors for the 
following 10^-10"^ years. This long term evolution of NS is independent of the cooling history during the 
first several hours when the interior cools down to temperatures of the order of 0.1 MeV. The thermal 
history of NS does strongly depend on the neutrino emission rates from dense matter during the neutrino 
radiation epoch t < 10^ — 10^ years. The neutrino emission rates in turn depend crucially on matter 
composition, elementary particle content, and condensed matter properties, such as superfluidity and 
superconductivity. Thus the studies of thermal evolution of neutron stars offer a unique tool to test 
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the physics of neutron star interiors. The measurements of the surface temperatures of young NS by 
X-ray satelhtes have the potential of constraining the properties of dense matter. The late-time t > 10^ 
thermal evolution of NS is dominated by cooling via the photoemission from the surface and heating due 
to conversion of rotational and magnetic energy into heat. This chapter reviews the neutrino radiation 
processes that are relevant for the neutrino radiation era. 

3.1 Classification of weak processes 

We start with a classification of the weak reactions and concentrate first on the nucleonic matter. It 
is useful to classify the processes by the number of baryons in the initial (or final) state. The simplest 
neutrino emission process that involves single baryon in the initial (final) state are 

n ^ p + e + u, p + e n + u, (157) 
n ^ n + u + i? (forbidden). (158) 

The first reaction is the charge current /3-decay (and its inverse). It is known in astrophysics as the Urea 
processes |187j . The Urea reaction is kinematically allowed in matter under /3-equilibrium if the proton 
fraction is sufficiently large |188l I189j . The threshold for the Urea process arises from the kinematical 
requirement of simultaneous conservation of momentum and energy in the reaction; a simple estimate 
for cold matter in /3-equilibrium shows that the proton fraction Yp > 11 — 14% for the Urea process to 

work pmimni . 

The second process - the neutral current neutrino pair bremsstrahlung is forbidden by the energy and 
momentum conservation. This statement is true if n refers to (quasi)particles whose spectral function 
is a delta function [cf. Eq. If one chooses to work with excitations that are characterized by finite 

widths, the reaction ()158|) is allowed |191j . The point is that the finite width incorporates multi-particle 
processes that we are going to include explicitly in the next to leading order of expansion. The process 
with two baryons in the initial (and final) states are the modified Urea and its inverse |192t I193j 



n + n—i'n+p + e + h', p + n^p + p + e + v, (159) 

n + p + e^n + n + v, p + p + e—^p + n + u. (160) 

and the modified bremsstrahlung processes 

N + N ^ N + N + iy + iy, Nen,p. (161) 



The modified processes are characterized by a spectator baryon that guarantees the energy and momen- 
tum conservation in the reaction. We can continue adding further spectator baryons on both sides of the 
reactions above, however the power radiated by higher order reactions drops dramatically for two rea- 
sons. First, if we consider the three-body processes, which are next in the hierarchy, the probability of 
scattering of three quasiparticles is suppressed compared to the two-body counterpart (Pauli principle). 
Second, adding an extra fermion in the initial and final state introduces a small factor T/Ep <^ 1 for 
each fermion, where Ep is the Fermi energy. Thus, going one step higher in the hierarchy suppresses the 
reaction rate by a small parameter {T / EpY- The relevant quantity for numerical simulations of neutron 
star cooling is the neutrino emissivity, which is defined as the power of energy radiated per unit volume. 
The emissivities of the processes above are ~ 10^"^ x T| for the reaction (jl57j) . Smod. p ~ 10^^ x Tg for 
the reactions ()159p and ()16Up . and ~ 10^^ x Tg for the reactions ()161|) : here Tg is the temperature 
in units 10^ K. 
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The general arguments above apply to the reactions in the hypernuclear matter under /3-equilibrium 
f Subsection I2.6|) . The charge current processes on hyperons (hyperon Urea processes) are 



n + e + p, ( ) ^ ( ) + e + ( ) ^ I S° I + ^ + ^ (162) 




and the reactions inverse to these. Note that some of the reaction involve violation of strangeness 
conservation by weak interactions with a change of strangeness by amount AS = 1. The neutrino 
emissivity of the hyperon Urea process has the same parametric dependence on the density of states 
and temperature as the nucleonic Urea process, however it is smaller than the nucleonic Urea process 
because the weak matrix elements are suppressed by factors of 0.01 - 0.6. Another important difference 
is that the threshold for the A° peu reaction, corresponding to the A° hyperon fraction Y\ ~ 0.0013, 
is much smaller than the threshold for the process ()157|) [194 . The remainder reactions have similar 
thresholds as the ordinary Urea and, being less effective, contribute a fraction of the total emissivity. 
The pair bremsstrahlung processes (|158|) on hyperons need to conserve the strangeness and are forbidden 
at the one-body level. The modified Urea processes involving hyperons can be written as 



Bi + B2^ B3 + B2 + e + u, (163) 

with the baryons Bi [i = 1,2,3) chosen consistent with charge conservation; the strangeness is either 
conserved or changed by amount AS = 1. The inverse of p63|) yields the same result as the direct 
reaction. The modified bremsstrahlung process can be written as . J95j 

5i + ^2 ^ ^3 + ^4 + z/ + z/, (164) 

where the weak interaction vertex, written symbolically a.s B ^ B + u + u, involves arbitrary octet 
baryons. (Note that contrary to ref. |195j . current work finds non- vanishing matrix elements for the 
A ^ A + u + u transition within the SU{Q) quark model for the baryon octet |196j ). The estimates 
of the reactions ()163p and ()164|) which are based on the one-pion-exchange (OPE) potential show that 
the emissivities of hyperonic processes are small compared to their nucleonic counterparts, contributing 
at most 2/3 of the total nucleonic emissivity. The OPE interaction adequately represents only the 
long-range part of the interaction. The short range repulsive components of the nuclear force reduces 
the rate of neutrino emission by modified processes by factors 4-5; since the reduction applies to the 
entire baryon octet the relative ratio of the nucleonic to hyperonic emissivities should not be affected. 
The exchange of strangeness carrying mesons, such as the i^'-meson, opens new channels for hyperonic 
modified processes and this will enhance the contribution of hyperons to the cooling rate. 

We mentioned above that the reaction (jl58p is forbidden for quasiparticles with 5-function type 
spectral functions ()45|) by the energy and momentum conservation. This constraint is lifted in the 
case where the baryons pair |197| 11981 1199j . The binding energy of Cooper pairs makes it energetically 
possible to create pairs in the inelastic processes of pair creation and annihilation: 



B{A) ^ B{A) + u + u, (165) 

where B{A) is a quasiparticle excitation of the superfluid state. The rate of the Cooper pair breaking 
and formation (CPBF) processes is of the order of (10^^ — 10^^) x Tg depending on the baryon pairing 
patterns |197l 11981 1199j . At extreme low temperature these processes are suppressed exponentially as 
exp(— A(0)/T), where A(0) is the zero-temperature pairing gap. 
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The neutral current one-body process ()165|) induced by the superfluidity have their charge current 
counterparts |2UUj . While the former vanish when the temperature approaches the critical temperature 
of superfluid phase transition, the emissivity of the latter process approaches the value of the corre- 
sponding Urea process. Put another way, the suppression of the Urea processes in the superfluid state 
is not restricted to the reduction of the phase space; the pair breaking can take place for these pro- 
cesses as well. In the extreme low temperature limit the Urea processes are suppressed exponentially as 
exp(— Amax(0)/T), where Amax(O) is the largest gap for fermions involved in the reaction. The modified 
processes are suppressed by factors exp{ — [Ab(0) + Ab'{0)]/T} where B and B' label the pair of the 
initial or final baryons. As in the case of the Urea process the largest gaps enter the suppression factor. 

If pions or kaons form a Bose-Einstein condensate (BEC) the meson decay reactions contribute via 
the reactions 201] 

n + (tt^) ^ n + e~ + Ue, (166) 
n + {K') n + + Ve- (167) 

The emissivities of these reactions are large compared to those of the baryonic processes above; for 
pions St, ~ lO^^Tg^ and for kaons Sk ~ lO^^Tg^; thus, a distinctive property of the models that feature a 
meson condensate is the rapid cooling. Note the kinematical differences in these reactions, since pions 
condense at finite momentum and in the P-wave, while kaons form a zero momentum condensate in 
the S-wave. 



3.2 Transport equations for neutrinos 

This section introduces the real-time formalism for neutrino transport \2Q2\ . We shall treat neutrinos as 
massless particles, since on the energy scales relevant for neutron star physics the masses of neutrinos 
are small. For massless neutrinos it is irrelevant whether neutrinos are Dirac or Majorana particles and 
their free particle Lagrangians are identical. The interaction Lagrangian for the charge current is 

= ^ E (^-"m^"'' + ^i^^^'') ' L-^ = ^lal,^La. Ll, = i'Lal,rLa, (168) 
* a 

where gw is charge current coupling constant, are the gauge vector boson fields, L^^ is the lepton 
current written in terms of left-handed chiral spinors of neutrino ip^j^ and lepton ip^^L^ ^ is the flavor 
index. The interaction Lagrangian for the charge neutral interaction is 

yE(^"^^'' + ^4^'')' = rLa7,rLa, = ^ll,{c. ' CAlVa, (169) 

a 

where gz = (7iy/cos 9]^, where 9^ = 28.7°, sin '^9w = 0.23. These coupling constants are related to 
the Fermi coupling constant by the relation Gp = {1 / \/i){gw /'^MwY = 1.166 x 10^^ GeV~^, where 
Mw — 82 GeV is the M^-boson mass (the Z-boson mass Mz = Mw/ cos 9w)- 

The theory of neutrino radiation can be conveniently formulated in terms of the real-time quantum 
neutrino transport, as discussed in Section |21 The neutrino Greens functions are written in the matrix 
form 

S%1,2) S<{1,2) \ _ [ (THm^)) -(V^(2)V^(1)> 



'^^^'^^-\s>ii,2) 3^(1,2) )-[ {^{i)m) {fmm} j' ^^^^^ 
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where are the neutrino field operators, = j^ip*, T is the chronological time ordering operator, and 
T is the anti-chronological time ordering operator; the indexes 1,2,... denote the space-time arguments. 
The neutrino matrix propagator obeys the Dyson equation 

5(1,2) = 5o(l, 2) + 5o(l, 3)0(3,2)5(2,1), (171) 

where 5o(l,2) is the free neutrino propagator and 5,^^(1, 2)5o(l, 2) = az5{l — 2), cr^ is the third 
component of the (vector) Pauli matrix, ft is the neutrino proper self-energy and we assume integration 
(summation) over the repeated variables. The self-energy f2 is a 2 x 2 matrix with elements defined on 
the contour. 

The set of the four Green's functions above can be supplemented by the retarded and advanced 
Green's functions which are defined, in analogy to (fTUI) and (fTT|) . as 

tS^{l,2) =e{h-h){{mM2)}), 2^^(1,2) = -^(t2-ti)({^(l),V^(2)}), (172) 

where { , } stands for an anti-commutator. The retarded and advanced Green's functions obey integral 
equations in the quasiclassical limit. By applying the Langreth-Wilkins rules fll5|) and ()16p to the 
Dyson equation !^ni\i we find the transport equation for the off-diagonal elements of the matrix Green's 
function 

[^3 - Re 3), 5>'<(3, 2)] - [Re S^{1, 3), fi>'<(3, 2)] 

= ^ {S>'<{1, 3), fi>'<(3, 2)} + i {fi>'<(l, 3), 5>'<(3, 2)} , (173) 

where [ , ] stands for commutator. In arriving at Eq. (jl73j) we assumed the existence of the Lehmann 
representation for the neutrino propagators; as a result we find Re = Re = Re S" and Re fl^ = 
Re = Re n. 

The neutrino dynamics can be treated semiclassically, by separating the slowly varying center-of- 
mass coordinates from the rapidly varying relative coordinates. Carrying out a Fourier transform with 
respect to the relative coordinates and keeping the first-order gradients in the slow variable we arrive 
at a quasiclassical neutrino transport equation 

i {ReS~^{q, x), S'>'<(g, x)}p^ + i {Re S{q, x), VL>^<{q, x)}pj^ 

= S>'<{q,x)n>'<{q,x) + n>'<{q,x)S>^<{q,x), (174) 

where q = {q, go) and x are the neutrino four momentum and the center-of-mass space-time coordinate, 
respectively, {. . .}p,b. is the four-dimensional Poisson bracket [cf. Eq. (jTTj) ]. To eliminate the second 
Poisson bracket on the 1. h. side of Eq. ()174|) we carry out a decomposition analogous to (j^Hj) with 
respect to the small neutrino damping: S^''^{q,x) = SQ'^{q,x) + Si'^{q,x), where SQ'^{q,x) is the 
leading (quasi-particle) and S^'^{q,x) is the next-to-leading order term. The quasiparticle part of the 
transport equation is then written as | l2U2j 

i{ReS-\q,x),S^'<iq,x)}p^=S>^<iq,x)n>^<iq,x) + n>^<iq,x)S>^<iq,x) (175) 

and describes the evolution of the distribution function (Wigner function) of on-shell excitations. The 
1. h. side of Eq. (|175|) corresponds to the drift term of the Boltzmann equation, while the r. h. side 
corresponds to the collision integral, where the self-energies Q^'^{q,x) are interpreted as the collision 
rates. The advantage of this form of the (generalized) collision integral is that it admits systematic 
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approximations in terms of the Feynman perturbation theory. The remainder part of the transport 
equation 

t{ReS-\q,x),S>^<{q,x)}p^ +t{Re S{q,x),Q>'<iq,x)}p^ = 0, (176) 

relates the off-mass-shell part of the neutrino propagator to the self-energies in a form of a local func- 
tional which depends on the local (anti-)neutrino particle distribution function and their coupling to 
matter. 



3.2.1 On-shell neutrino approximation 



The on-mass-shell neutrino propagator is related to the single-time distribution functions (Wigner 
functions) of neutrinos and anti-neutrinos, fu{q,x) and fp{q,x), via the ansatz 



Soiq,x) 



ITT 



S (go - ujuiq)) fu{q, x) - S (go + uj^iq)) (1 - fu{-q, x)) 



(177) 



where u!i^{q) = \q\ is the on-mass-shell neutrino/anti-neutrino energy. Note that the ansatz includes si- 
multaneously the neutrino particle states and anti-neutrino hole states, which propagate in, say, positive 
time direction. Similarly, the on-shell propagator 



S>{q,x) 



in 4 



5 (go - ujy{qy) (1 - /^(g, a;)) - 5 (go + ujy{qy) fp{-q, x) 



:i78) 



corresponds to the states propagating in the reversed time direction and, hence, includes the anti- 
neutrino particle states and neutrino hole states. 

To recover the Boltzmann drift term, we take the trace on both sides of the transport equation ()174|) 
and integrate over the (anti-)neutrino energy go. The single time Boltzmann equation (hereafter BE) 
for neutrinos is obtained after integrating over the positive energy range: 



dt + dgU^{q)d^ fu{q,x) 



dqo 
2tt 



Tr [n< (g, x)S> (g, x) - n> (g, x)S< (g, x)] 



(179) 



a similar equation follows for the anti-neutrinos if one integrates in Eq. p74j] over the range [—oo, 0]. 

The different energy integration limits select from the r. h. side of the transport equations the 
processes leading to modifications of the distribution functions of (anti-)neutrinos. The separation 
of the transport equation into neutrino and anti-neutrino parts is arbitrary. It is motivated by the 
observation that the fundamental quantities of neutrino radiative transport, as the energy densities or 
neutrino fluxes, can be obtained by taking the appropriate moments of BEs and these quantities are 
not symmetric with respect to the neutrino/anti- neutrino populations in general. 



3.2.2 Collision integrals 

The diagrams contributing to the neutrino emission rates can be arranged in a perturbation expansion 
with respect to the weak interaction. The lowest order in the weak interaction Feynman diagrams which 
contribute to scattering, emission, and absorption processes are shown in the Fig. EH The corresponding 
transport self-energies are read-off from the diagram 

-^^]>'<(gl,x) = y"^|^(2vr)V(gi-g2-gX,^^o<fe,x)^riXA<(g,x), (180) 
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Figure 20: The neutrino self-energies for charge and neutral current processes. The dashed 
and double-dashed curves correspond to the free and full neutrino propagators (reverting the 
time-direction one finds the propagators for anti- neutrinos). The solid line is the electron 
propagator. The loop is the baryon polarization tensor involving baryons and B^- The 
wavy lines correspond to the and Z° boson propagators. The incoming and outgoing v 
propagators are shown for clarity and are not included in the self-energy. 



where nj^j^^(g) are the off-diagonal elements of the matrix of the baryon polarization tensor, F^^ is the 
weak interaction vertex. The contact interaction can be used for the energy-momentum transfers much 
smaller than the vector boson mass, q <^ M2,Miy, in which case the gauge boson propagators are 
approximated as 

'^^'^ = a^-M^ ~M^' ^^^^^ 

where Qf^u is the metric tensor. Let us first concentrate on the BE for neutrinos. Define the loss and 
gain terms of the collision integral as: 

/>'<(q,x) = r^Tr [n>'<iq,x)S>'<{q,x)] . (182) 

Jo ^TT 

Substituting the self-energies and the propagators in the collision integrals we find for, e.g., the gain 
part: 

luiQux) = -i [ ^Tr| /" ^ ^ f {27ry S^{qi - q2 - g)F^ ^ \s (go2 - ^,.(^2)) 1^2, x) 

- S (go2 + ^^.(92)) (1 - fp{-q2,x)) — {qio - uj^{qi)) (1 - /^(gi, x)) (g, x) \. 

J uj^[qi) ) 

(183) 

The loss term is obtained by replacing in Eq. ()183|) the neutrino Wigner functions by the neutrino-hole 
functions fu{q,x) ^ (1 — fu{q,x)) and the anti-neutrino-hole Wigner functions by the anti-neutrino 
functions (1 — fu{—q, x)) — > /p(g, x). The terms proportional (1 — fu)fu and (1 — /^)(1 — /p) in the gain 
part of the collision integral, I^{q), correspond to the neutrino scattering-in and emission contributions, 
respectively. The terms proportional fu{l — fu) and f^fp in the loss part of the collision integral, I{^{q), 
are the neutrino scattering-out and absorption contributions. 

The loss and gain collision integrals for the anti-neutrinos can be defined in a manner, similar to 
the case of neutrinos, with the energy integration spanning the negative energy range 

/p>'<(q,x) = I' ^Tr[^]>'<(g,x)So>'<(g,x)]. (184) 
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Using the above expressions for the self-energy and the propagators, we find, e.g., for the gain term: 



° dqio 
^ 2tt 



Tr 



oo ^4 



5 (go2 - ^^u{q2)) fu{q2,x) 



- 6 (go2 + ^^u{q2)) (1 - /p(-g2, x)) 



rt'^^5(gio + ^.(qi))/p(-gi,x)n> (g,a;) (185) 



The loss term is obtained by making replacements in Eq. ()185j) analogous to those applied to Eq. 
(1 18 Hp . The terms proportional f^fu and /p(l — /p) in the gain part of the collision integral, Ip{q), 
then correspond to the neutrino absorption and scattering-out contributions. The terms proportional 
{l—fp){l—f^) and {l—fu)fu in the loss part of the collision integral, Ip{q), are the neutrino emission and 
scattering- in contributions, respectively. Note that, when the neutrinos are in thermal equilibrium with 
the baryons, the collision integrals for the scattering-in/scattering-out and for the absorption/emission 
cancel. Under the conditions of detailed balance the (anti-)neutrino distribution function reduces to 
the Fermi-Dirac form. 



3.2.3 Neutral current processes (bremsstrahlung) 

The neutrino-pair emissivity (the power of the energy radiated per volume unit) is obtained by mul- 
tiplying the left- hand- sides of the neutrino and anti-neutrino by their energies, respectively, summing 
the BEs, and integrating over a phase space element: 



JtJwf ^^^^^^ ^ ^'^^^^ '''^^^^ ^Iw)~' ^^'^''"^^^ ~ ^^'"'^^^^ ^^^^^ 



where in the collision integrals we kept only the terms which correspond to the processes with the 
neutrino and anti-neutrino in the final state (bremsstrahlung) 

{2iry J {2TTy2uj^{qi) {27Ty2uj^{q2) {2tc)'^ 

+ ^M2) - go)^.(gi) [1 - fuMqi))] [1 - hMq2))] A^\qu g2)n;:,<(g, x),(187) 

and A^'^ = Tr [7^(1 — 7^) 4il'^{^ ~ 7^) ^2]- The collision integrals for neutrinos and anti-neutrinos can 
be combined if one uses the identities Ii'^-^{q) = n^^(— g) = 2igB{qo)lio^^^xil)'^ here gsilo) is the Bose 
distribution function and 11^;^ (g) is the retarded component of the polarization tensor. With these 
modifications the neutrino-pair bremsstrahlung emissivity becomes |19H 11981 l2U2j 

G f d^q2 f dSi f A 



2V2J yj {2Txf2uo,{q2) J (27r)32a;,(gi) i (27r)^ 



(27r) 5\qi + q2- q)6{uj^{qi) + ^^(^2) - go) K(qi) + ^uiq2 



9B{qo) [1 - fuMqi))] [1 - /p(cu,(q2))] A'^'(gi,g2)Imn^,(g). (1^ 



The symbol Im refers to the imaginary part of the polarization tensor's resolvent and the /-sum is over 
the three neutrino flavors. We note that Eq. ()188j) is applicable for arbitrary deviation from equilibrium. 
Therefore Eq. ()188|) is applicable beyond the boundaries of the linear response theory or the S'-matrix 
theory which explicitly resort to the equilibrium properties of the system as a reference point. 



54 



3.2.4 Charged current processes ((3-decay) 

In the case of charge current processes there is a single neutrino or anti-neutrino in the initial and 
final states. It is sufficient to compute, say, the direct /3-decay and multiply the result by a factor 2 to 
account for the inverse process. The anti-neutrino emissivity is obtained from the counterpart of Eq. 
flT7n|) for anti-neutrinos PHH] : 

In full analogy to the charge neutral current interactions we obtain 

f gV f d^qi r d^q2 /■ ,4 

S(u,, + w„ - q„) u,,{q2)9B{qo) |1 " M'^.)] |1 " Uf )] A"*(9i, ?2) Im n« (g), (190) 

where Ue is the electron energy, G = Gpcos 9c, and 9c is the Cabibbo angle (cos 9c = 0.973). Note 
that in cold matter (temperatures T < 5 MeV) neutrinos propagate without interactions and to a good 
approximation f,^, fp <^ 1. The properties of matter to which neutrinos couple are encoded in the 
polarization tensors according to Eqs. ()188p and ()190p . 

3.3 Polarization tensors of hadronic matter 

The neutrino emission rates depend on the response of hadronic matter to the weak probes in the 
time-like domain. The polarization tensors appearing in Eqs. pSSj) and ()190|) can be viewed as self- 
energies of and bosons. The neutrino self-energies in Fig. 1201 are then re-interpreted as the Fock 
contributions due to the exchange of renormalized gauge bosons. The picture adopted in Subsec. 13.2.21 
interprets the same diagrams as second order in weak interaction Born self-energies of neutrinos. Both 
interpretations are equivalent of course. 

The classification of the reactions by the number of baryons participating in the reaction (Subsec. 13. 1|) 
translates into expansion of the polarization tensor in particle-hole loops. The one-body processes ()157j) 
and p58j) are described by the one-loop polarization tensor, the two-baryon processes, e.g. Eqs. ()159j) 
and ()161|) are described by the two-loop polarization tensor, etc. This is the case if one works with 
well-defined quasiparticles described by the 5-function spectral functions. For general forms of spectral 
functions with finite width the situation is more complex: the loop expansion in the particle-hole channel 
can still be carried out, however the width of the spectral function should not contain resummation in 
this channel to avoid double counting. 

3.3.1 One-loop processes 

The one-loop polarization tensor is shown in Fig. |^ the corresponding analytical expression is 

^n<.(g) = J ^^0^jT[T,G<{p)TM^{p')]{2rrmq + p' -p) =T,X{ql (191) 

where the charged current weak interaction vertices are F^ = Gp'^fii^ — QaI^) with qa = 1.26 being 
the axial coupling constant, the tensor T^^, = —G^gl where gx = ^ ior fi = u = and gx = gA for 
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Figure 21: The one- loop polarization tensor for charged current process. The wavy lines 
correspond to the propagators, the solid line to the baryonic propagators. 



jj, = u = 1,2,3; the loop integral can be computed analytically in the quasiparticle limit 

1 + exp - /ib)] 



1 + exp [-(3 (e - I^b) - M 



"B' 



2nP\q 



-L{u,q), (192) 



where is the effective mass of a quasiparticle of baryon type B, ^ = /2m — andp = {m*^/q){uj — 
jJ'B + I^B' — 9^/2^^), where we assume that ~ m|j,. This result can be substituted in Eq. ()190|) to 
obtain, for example, the emissivity of the Urea process n — > p + e + z/ |2(JUj 

ep = eo y dy g{y) L{y,pFe) j dzzfe{z-y), eo = (1 + 3^^) ^^^^^ , (193) 

where ppe and /e are the Fermi-momentum and distribution function of electrons and y = (3uj. In the 
zero temperature limit L{y,pFe) = yG{.—K)^ the integrals in Eq. ()193|1 can be performed analytically 
and one recovers the zero-temperature result of Lattimer et al. |189j . The zero temperature 6'-function 
can be rewritten as 0{ppe + PFp — PFn) |189j which tells us that the "triangle inequality" Pfs + Pfp > Pfu 
must be obeyed by the Fermi-momenta of the particles for the Urea process to operate. 

The emissivity of the Urea processes p93p scales as since the process is second order in the 
weak interactions; its linear dependence on the effective masses of participating baryons arises through 
their density of states u ~ m*pF/'n'^, where pp stands for the Fermi momentum; the density of states of 
massless electrons contributes the factor and a factor 1/pFe arises from the momentum conservation 
condition. The temperature dependence can be understood from the dimensional analysis of the reaction 
rate and arises as follows |19Uj : each degenerate fermion being confined to a narrow band ~ T around the 
Fermi surface contributes a factor T, the final state anti-neutrino contributes a factor T^; an additional 
factor T is due to the fact that we are interested in the energy rate and a compensating factor of 1/T 
arises due to the energy conservation constraint. 

Eq. fll93|) can be adapted to describe the Urea processes ()163|) which involve hyperons |194j . Among 
these the process A'' —>■ p+e~ + u is potentially important. The corresponding triangle inequality follows 
from the theta function 9{ppf, + ppp — Pfa)- Since typically ppe ~ Pfp Pfu a small fraction of A's 
is sufficient for the reaction to operate. Although it is less effective than the nucleonic Urea process, 
since it involves a change of strangeness and its matrix element is proportional to sin "^Oc, it is still as 
efficient as other exotic cooling channels e.g. the pion decay p66p or kaon decay p67|) . 

For identical baryons the one-body bremsstrahlung process ()158|) vanishes, as can be seen from 
Eq. (|192j) in the limit of equal chemical potentials (note that in the non-relativistic kinematics spurious 
terms remain, which vanish exactly if the proper relativistic kinematics is used). However there are 
cases where the one-loop bremsstrahlung is possible because the baryons are embedded in a mean-field; 
an example are the CPBF processes which arise due to the pairing mean field (see subsection 13. 5p . Yet 
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another possibility arises when the baryons are coupled to external gauge fields. Canonical neutron 
stars support magnetic fields, of the order of 10^^ — 10^^ G. A separate class of neutron stars, know as 
magnetars, are believed to support fields that are much larger, of the order of i? ~ 10^^ — 10^^ G |2U3j . 
The interaction of the baryon magnetic moment with the applied field induces a splitting in the energy 
of spin-up and spin-down baryons of the order of fx^^^B, where B is the magnetic induction, fi^^^ is 
the fermion magnetic moment. For neutrons /i'-"'' = gnfJ-N, where the gyromagnetic ratio gn = —1.913; 
for protons fi^^^ = gp^N, gp = 2.793, where /xat = eh/2mp = 3.152 x 10~^^ G~^ MeV is the nuclear 
magneton. Thus, we can expect non-vanishing neutrino bremsstrahlung as a result of the Pauli spin- 
paramagnetic splitting whenever the 2fi^^^B ~ T |2()4j . For neutrons the Pauli paramagnetism is the 
only effect that affects the quasiparticle spectrum, which is written as (to the leading order in B'^/rriN) 

Bnis) = {ml + p2 - 2sm„/x(")5)l/^ (194) 

where the spin projections on the magnetic axis are s = ±1. The quasiparticle spectrum of charged 
particles includes in addition to the Zeeman splitting the Landau quantization of orbits in the plane 
transverse to the direction of the applied field 

= [^^ +pl + i2N + l- sgp)eBY/^, (195) 

where N labels the Landau levels. The emissivity due to the n ^ n + u + u process is obtained as |204j 

2(27r)5 Jo ^ey-ljQ \ e-l^^—v + l e^PE+-y + lJ' ^ ' 

where E± = (rrij^ + p| — sgneBY^"^ and p± are the upper and lower bounds on the momenta imposed 
by the integration over the angle between the particle momentum and the momentum transfer in the 
process. The emissivity involves only axial currents because the process requires a spin-fiip whereby 
a quasiparticle is transfered from one Fermi surface (of the spin-up population) to the other surface 
(spin-down population). The dependence on the effective masses arises from the phase space integration 
which introduces a density of states per particle. The analysis of the temperature dependence is similar 
to the case of the Urea process, with the only difference that the transfered momentum is of the order 
of T (rather than ppe) therefore an additional factor T appears. Note that this bremsstrahlung process 
vanishes in the zero field limit. Several authors considered the Urea process in strong magnetic fields, 
where the effects of Landau quantization change the qualitative picture by removing the "triangle 
inequality" constraint [M IM M • 

The temperature region where the one-body neutrino-pair bremsstrahlung is important increases 
with increasing magnetic field (see Fig. [221) • The pair bremsstrahlung from neutrons is efficient whenever 
\^n\B ~ T, since then the energy involved in the spin-fiip is of the same order of magnitude as the 
thermal smearing of the Fermi surface. The temperature at which neutrino-pair bremsstrahlung from 
neutrons becomes comparable to the competing processes roughly coincides with this condition. For 
lower temperatures the emissivity drops exponentially, because the energy transfer becomes larger 
than the thermal smearing. Neutrino-pair bremsstrahlung from protons is important when iipB ~ T. 
The emissivity due to the protons increases faster than the emissivity due to the neutrons with the 
temperature, since the smearing of the proton transverse momenta provides an additional relaxation on 
the kinematical constrains. According to Fig. |221 the emissivity of the modified Urea process is larger 
than that of neutrino-pair bremsstrahlung from neutrons and protons at high temperatures, mainly due 
to the different temperature dependencies of these processes (oc T'' for the one-body bremsstrahlung as 
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Figure 22: The emissivities of various processes vs the temperature for fields B = 10^^ {left 
panel) and B = 10^"^ G {right panel). 



compared to oc for the modified Urea). However, for temperatures that are smaller than fipB ~ T 
the emissivity drops just as for neutrons - exponentially. In the case of a superstrong magnetic field 
B > 10^^ G the large uncertainty in the transverse momenta of the protons and electrons allows the 
direct Urea process to occur, and its emissivity dominates the emissivity of any other process. 

3.3.2 Two-loop processes 

The nuclear interaction enters the quasiparticle loop expansion at the second order. To compute the 
emissivity we need a model of nuclear scattering in background medium. The form of the nuclear 
interaction depends, of course, on the nuclear matter model one works with. Below we will give a specific 
example of the computation of the neutral current bremsstrahlung process n + n^n + n + u + u |2(J2j . 
The effective particle-hole interaction can be represented by pion-exchange at long distances and contact 
Fermi-liquid interaction at short distances |207j 

V[ph]{k) = (^] (o-i ■ k) D\k) {a2 ■k) + fo + h{(T, ■ (T2), (197) 

where f-,, is the pion decay constant, m^r is the pion mass, D^{k) is the on-shell causal pion propagator, 
/o and /i are the coupling parameters of the Fermi-liquid theory, cr is the vector of the Pauli matrices. 
The form of the ph interaction is suitable when the scales in the problem can be separated with respect 
to the Compton length of the pion A^^ = = 1.4 fm. If the system is characterized by scales \ 
(e.g. is sufficiently dilute) the only relevant dynamical degree of freedom is pion and the rest of the 
nuclear interaction can be approximated by constants. To obtain values of neutrino emissivities that are 
consistent with those computed from the nuclear T-matrix the p meson exchange needs to be included 
exphcitly in Eq. (IT!?7j) . 
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Figure 23: The two-loop polarization tensor. The wavy lines correspond to the propa- 
gators, the solid lines to the baryonic propagators, the dashed lines to nuclear interactions. 



The topologically non-equivalent two-loop diagrams of our theory are shown in Fig. 1221 The ana- 
lytical expression, corresponding to the first (from left to right) diagram in Fig. EHl is 



i=l 



dk 

(27r)«5(g + - k - p^)5{k + p2 - pi)Tr [Vik)G<ipi)V{k)G>{p2)] 



(27r 



Tr [r^G^(g + p,)Vik)D%k)G<ips)Vik)D''ik)G''iq + p,)T,G>ip,)\ , (198) 

where V{k) is the strong interaction vertex, which can be read-off from Eq. ()197|1 . The contribution 
of this diagram is readily recognized as a propagator dressing in the ph channel by a self-energy which 
involves an excitation of particle- hole collective mode. The second diagram in Fig. 1221 corresponds to a 
vertex correction in the ph channel by an effective interaction, which incorporates an intermediate the 
particle- hole collective mode excitation. The third diagram in Fig. 1221 may be interpreted as a particle- 
hole fluctuation. These diagrams for model interaction ()197|1 show the following features: (i) the vector 
current contributions from the first and second diagrams mutually cancel; (ii) the third diagram does 
not contribute because the axial-vector contribution involves traces over odd number of cr-matrices 
and the vector-current contribution is canceled by an equal and of opposite sign contribution from the 
diagram which is generated from the third diagram by flipping one of the loops upside-down; (iii) all 
contributions due to the Fermi-liquid interaction cancel after summing the first two diagrams. Note 
that the exchange diagrams are generated from the direct ones by means of interchanging the outgoing 
propagators in a strong vertex. 

The causal propagator in Eq. ()198|) have the following general equilibrium form 



G^ip) = "-^V^^ tanh ( ^) , (199) 

where tanh {uj/2) = [1 — 2/^^(0;)] and = p^/2m-|-Re S(p). The arguments of the causal propagators in 
the diagrams in Fig. (221 contain the external four momentum q and the propagation between the strong 
and weak vertex it describes is off the mass-shell. The off-shell dependence of the propagator can be 
simplified by expanding with respect to small parameter f ^ 1 which is the characteristic velocity of a 
baryon (the velocity of light c = 1 and we use non-relativistic kinematics) 

d 

{uj + Ep) - ~ - p ■ q/m - q— Re S(p) - ~ u;, (200) 

to the leading order. For off-shell energies not far from the Fermi energy the quasiparticle damping 
is an even function of the frequency 7(— cu) = 'j{uj) (this observation is exact in the phenomenological 
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Fermi-liquid theory and is confirmed by microscopic calculations). Then 



G%±u:,p) = ±-, T tanh V (201) 

Note that the acausal propagator is obtained from the above through complex conjugation [G'^{p)]* = 
—C^ip). We conclude that the propagator ()201|) is odd under the exchange of the sign of tu, a property 
which is important for the vector current conservation. We now write down the neutrino emissivity for 
the n + n^n + n + v + v process which is the sum of diagrams in Fig. |2S1 

GWa f ^ V f — V / = 5.5 X 1019 h Tl (erg cm'^ s'^) (202) 



5(27r)9^^^V"^^/ V"^ 
where T9 is the temperature in units of 10^ K, is the integral / in units 10^ defined as 

noo poo poo 

1= / dyy'\G%y)\''Q{y) dxx'\D'^{x)\^ dz g{z) g{y - z)L{z,x)L{y - z,x). (203) 

Jo Jo J-oo 

The temperature dependence in Eq. (I2()2|l can be understood by comparing the modified processes 
to their one-body counterparts discussed above. The additional fermion appearing in the initial and 
final state introduces an additional factor (T/Ep)^. Since this argument applies equally well to the 
modified Urea process, we conclude that its emissivity scales with temperature as T^. The scaling of 
emissivity with the effective mass arises due to its dependence on the density of states of the initial and 
final baryons. Each strong interaction vertex introduces a power of the pion-nucleon coupling, therefore 
ej,p oc (see Fig. 1221) • The simplest approximation to the pion propagator in Eq. ()203|) is to neglect the 
pion self-energy and approximate it as D'^{k) = [k"^ + m^]"^. The free-space approximation should be 
valid in the vicinity of the nuclear saturation density and below. The softening of the one-pion exchange 
(a precursor of the pion-condensation) increases the neutrino emissivity by large factors, see for details 
ref. |2()lj . The Pauli blocking of the final state neutrino at finite temperatures is taken into account 
by the function Q{y); In the dilute (anti-) neutrino limit Pfi^f ^ 1 (where fi^^ is the chemical potential 
of neutrinos of flavor /) Q{y) = 1. This is the case below the temperatures where the neutrinos are 
trapped. In the low-temperature limit C{z) = z and the ^-integration decouples from the x-integration. 
Letting 7(0;) — »• (quasiparticle limit) one finds that G'^{uj) = uj~'^. Then, the z integration can be 
carried out analytically upon dropping the wave-function renormalization contribution: 

dz g{z) g{y - z) z {y - z) = ^ J > . (204) 

6 {ey - 1) 

After these manipulations Eq. ()2()2j) reduces to the quasiparticle result of ref. |2()7j . It should be 
noted here that the OPE approximation to the nucleon-nucleon amplitude is not justified in dense 
matter from the numerical point of view, and one should include other mesons to take into account 
for the intermediate range attraction and short-range repulsion. In particular the inclusion of the p 
meson repulsion modifies the meson propagator to I208lj D'^{k) = [fc^ -|- rri^]^^ — Cp[k'^ + ^p]"^, where 
Cp = 1.67 and rup ~ 600 MeV. Such a correction substantially improves upon the OPE result and a 
quantitative agreement is achieved between the vrp-exchange model and full T-matrix calculations 
or low-momentum interactions 
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Figure 24: The on mass-shell damping 7 [upper panel) and the real part of the self-energy 
[lower panel) as a function of the particle momentum for neutron matter at density po = 0.16 
fm~^. 



3.3.3 Landau-Pomeranchuck-Migdal effect 

In this subsection we explore the effect of the finite width in the causal propagator ()20H) which de- 
scribes the off-shell propagation between the weak and strong vertices. The off-shell propagations is 
characterized by a length (or time) scale known as the formation length (time) first introduced in the 
context of bremsstrahlung of charge particles passing through matter by Ter-Mikaelian |2111 1212j . In 
our context the gauge boson energy is soft uj <^ Ep and can be associated with the formation length 

h 

If = —vp = TfVp, (205) 
to 

where vp is the baryon Fermi- velocity. The formation length is the distance that a particle covers during 
the emission of the gauge boson; if vp is large and uj small, // can be very long. This observation was the 
basis for the suppression calculation by Landau and Pomeranchuk |213j and Migdal |214j in the context 
of high-energy electrons radiating photons. There is another scale in the problem - the mean-free-path 
of a quasiparticle, 

Imfp = TynfpVp = [na)'^ ("Y") 

where a ~ 40 mb is the baryon-baryon cross-section (mb = 10~^^ cm^), n is the baryon number density. 
In the limit l^fp ^ If the radiation is from a well defined source - the environment has no influence 
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on the radiation. In the opposite hmit Imjp ^ the baryon-baryon interaction need to be included in 
the renormahzed propagators, which are now the fundamental degrees of freedom to treat the radiation 
process. In the case where Z^/p ~ the radiation is suppressed since the baryon-baryon collisions 
interfere (destructively) with the emission process. 

In our context the Landau- Pomeranchuk-Migdal (LPM) effect is described by the finite width of the 
quasiparticles in Eq. ()2UH) . The finite width not only suppresses the radiation, but also regularizes the 
infrared divergence of the radiation cross-section; in the case of the neutrino-bremsstrahlung processes 
this divergences is absent even in the zero width limit because of high powers of uo appearing in the 
numerator of the emissivity. 

The width of the quasiparticle propagators can be parameterized in terms of the reciprocal of the 
quasiparticle life time in the Fermi-liquid theory: 



where a is a density dependent phenomenological parameter. Fig. |^ illustrates the real part of the 
self-energy and the damping in neutron matter computed within the finite temperature version of the 
BBG theory. In this version of the theory the K matrix is complex valued |2U2j . 
The emergent neutrino spectrum can be characterized by the spectral function 



The dependence of the integral on the (dimensionless) neutrino frequency y = jSu at T = 20 MeV and 
the saturation density po is shown in Fig. [JHlfor three cases: the limit of vanishing width {dashed line), 
including the leading order in 7 contribution in the width [i. e. the first term in Eq. ()2()l|l ] {dashed- dotted 
line) and the full non-perturbative result {solid line). The energy carried by neutrinos is of order of 
uj ~ 6T in all three cases since the average energy carried by each neutrino (non-interacting relativistic 
massless particle) is 3T. The finite width of propagators leads to a suppression of the bremsstrahlung 
rate. Keeping the full non-perturbative expression for the causal propagators enhances the value of the 
integral with respect to the leading order perturbative result. The LPM effect sets in when w ~ 7. 
As neutrinos are produced thermally, the onset temperature of the LPM effect is of the order of 7. 
Equation ()207|) shows that the value of the parameter a controls the onset temperature which turns out 
to be of the order of 5 MeV. 

An alternative to the two-loop calculation outlined above is the computation based on a one-loop 
polarization tensor with fully off-shell propagators, as suggested in the case of photoemission in ref. |215j 
and neutrino emission in ref. |191j . However, dressed propagators need to be supplemented with dress 
vertex functions that satisfy the Ward identities in self-consistent manner; the problem of dressed 
vertices is discussed in ref. |216j . Multi-pair excitations processes, that are relevant to the description 
of the LPM effect within the Landau- Fermi liquid theory are discussed in Refs. |2171 1218j . 

3.3.4 Soft-neutrino approximation 

The soft neutrino approximation arises in the models that are based on either the free-space or medium 
modified T-matrix theories |219| 122m I221j . We have seen that the typical energy carried by neutrinos 
is of the order of several MeV, which is small on the nuclear energy scales of tens of MeV. The neutrino 
bremsstrahlung process is therefore called "soft" . Note that this is not the case for the Urea processes, 
where the electron energy is of the same order of magnitude as the neutron energy. In the case of the 
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Figure 25: The neutrino spectral function ()208p at the temperature T = 20 MeV and density 
Po = 0.16 fm~'^. The dashed curve is the zero width hmit, the dashed-dotted curve includes 
only the leading order in 7 contribution from the causal propagator, the solid curve is the 
full non-perturbative result. 



bremsstrahlung one can apply the ideas underlying the Low theorem for the photon bremsstrahlung 
which states that the bremsstrahlung amplitude to the leading 0(x~^) and next-to-leading 0(1) order 
in the expansion with respect to x = ^/Ej, where uj and Ep are the characteristic energies of neutrinos 
and baryons, is determined by the non-radiative cross- sect ion. The weak matrix element is written as 



M''^ = T{p,p',;p, - q,P2)G{p, - q)T^ + ^^G{p[ + q)T'{p[ + q,p',;puP2) + (1 2), (209) 



where T is the scattering T-matrix, G{p) = mA+/(p • q) is the free-space Green's function with the 
positive energy projector defined as A+(p) = (j^ + m)/(2m), = [Gp / \/^)'yv{cv — CA){Ta/2) is the 
weak interaction vertex and t"- is the isospin operator. Next, expand the T- matrices in ()209|) around 
their on-shell matrix Tq 



The expansion of Eq. ()209|1 can be divided into the external and internal contributions, depending 
whether the neutrino is emitted from the external leg of the T-matrix or from the internal interaction 
line (see Fig. I26|). For the vector current the external term can be computed explicitly while the form 
of the internal term is fixed from the requirement of the vector current conservation q'^Ai^'"' = 0. For 
the axial vector case one needs to take into account the concept of partially conserved axial currents 
(PCAC) which leads to the condition |2ini 
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Figure 26: Illustration of the external {left graph) and internal {right graph) contributions 
to the neutrino bremsstrahlung. The shaded block corresponds to the T-matrix, the dashed 
line to the weak interaction, the arrows to the participating particles. 



where Ai'^ is the pion emission matrix element. To the leading order in the expansion with respect to 
the small parameter x, i. e. 0{x~^) the vector and axial vector matrix elements read 



GpCy I pi^ ^ ^ 



^V,a ^ I _y^:^^a ^ ^a^y ) + {1 ^ 2}, (212) 

" 2V2 \ Pi-q p'l-q J 

. rnGpCy [ ^ A+(pi) ^, , ^, ^, ^^ K+{p{) ^\ , _ ^1 

^ V2 V Pi ■ ? Vi-q J 

It can be verified that the vector current is conserved. For matrix elements with the nucleons on the 
mass shell, the most general nonrelativistic charge-independent nucleon-nucleon amplitude contains 
central, spin-spin, tensor, spin-orbit, and quadratic spin-orbit terms. For fixed isospin it can be written 

T = Tc + Ts(Ti ■ 0-2 + TtSu + TsoL ■ S + TqQ^ (214) 

where the tensor and quadratic spin-orbit operators are defined as = Scri ■ r<j2 ■ r — cr^ ■ cr^ and 

Q\2 = (CTI ■ L(72 ■ L + 0-2- L(Ti ■ L)/2. 

Consider the process n + n-^n + n + i' + u. In the non-relativistic limit p ■ q = p ■ q — fnuj ~ mu 
to leading order in Vp/c and the vector current contribution ()212j) vanishes. The contribution from the 
axial-vector current is given by the commutator 

-Mo=0, M^ = ^^[T'^^5], (215) 

where S = (cri + a2)/2, is the spatial component of the vector A4^, T"'" is the neutron-neutron 
scattering T-matrix and the axial coupling constant qa = 1.26. The commutator ()215p is non-zero 
only for the tensor, spin-orbit and quadratic spin-orbit terms in the expansion ()214|) . therefore only 
spin-triplet nucleon-nucleon partial waves contribute. The dominant tensor force is much stronger in 
the np system and despite the fact that the proton fraction is small the process n + p^n + p + u + i' 
can gain significance. The vector and axial- vector matrix elements of this process in the non-relativistic 
and soft-neutrino approximations are 

M^ = -^.M'', M"" = -%^^—^-T-^, (216) 
LU 2^2 Lu m 

= 0, M^ = ^^[T^p^ S-], (217) 
v2 "-^ 
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where k = p — p', S~ = (cri — cr2)/2, Cy = —1, Cy = 1 — iSin^O^r. The soft neutrino approximation is 
not vahd for the modified Urea process (jl59|) as the energy transfer in the reaction could be large, of 
the order of the neutron Fermi-energy. The neutrino emissivities evaluated with the matrix elements 
quoted above (which thus include the full T-matrix) are reduced by factor 4-5 compared to the results 
obtained with the one-pion-exchange amplitudes |219| I22U[ I22H I^U^ . but are of the same magnitude if 
the OPE is supplemented by repulsive p-meson exchange. 



3.4 Graviton emission in Kaluza-Klein theories 

We have seen that the charged current two-loop process - the modified Urea reaction - is dominant 
among these type of reactions. There is a distinct physical situation where only the charge-neutral 
current processes contribute - the case of gravitational bremsstrahlung f 2211 1221^ . This process arises 
within the models where the standard model fields live on a four-dimensional manifold, while gravity is 
allowed to propagate in {n + 4) dimensions, where n is the number of extra dimensions. Such schemes 
provide a natural explanation of the hierarchy problem in particle physics; the gravity appears to be 
much weaker that the remaining three fundamental forces because it is diluted by its extension in extra 
dimensions. This picture fits into the traditional Kaluza-Klein (KK) theories which contain the usual 
four dimensional space-time manifold plus additional compact dimensions which form an unobservable, 
small manifold (which until recently was believed to be Planck size ~ Mp\ where Mp = 1.2 x 10^^ 
GeV). Recent models of extra dimensional gravity propose that there are n extra compact dimensions, 
all of which are about the same size R and R is much larger than the Planck scale, possibly as large as 
a millimeter |22H| I224j . The size R of the extra dimensions is given in terms of the Planck mass Mp 
and an effective mass M, which is taken of the order of 1 TeV, as 

^" ~ TT^- (218) 

At scales of the order of R the Newtonian gravity is expected to fail. If n = 1 then for M ~ 1 TeV 
one finds R ~ 10^° km, which implies that there must be deviations from Newtonian gravity over solar 
system distances. However, if n = 2 then R < 1 mm. Since gravity has not been tested at distances 
smaller than millimeter, large extra dimensions are consistent with present experimental knowledge. 
An interesting consequence of this theory is that a 4 + n dimensional graviton can propagate in extra 
dimensional space, while the standard model particles are confined to the four dimensional space. 

The size of the extra dimensions can be constrained if there is an evidence of missing energy in 
astrophysical processes such as supernova explosions. Since the models of supernovae based on the 
standard physics explain the duration and energy of the neutrino pulse observed from SN1987A, any 
mechanism that drains sufficient energy from the core of the supernova will destroy the agreement. The 
bremsstrahlung of gravitons in the nucleon-nucleon collision was suggested as such a mechanism |223j . 
The interaction between gravitons and dilatons and the standard model particles is described by the 
Lagrangian density 

1/2 



(219) 



where h^^'^ and 0^ are the graviton and dilaton fields, j is an n-dimensional vector representing the 
momentum of the mode in the extra dimensions. The differential rate at which the KK particles escape 
into extra dimensions can be related to the on-shell T-matrix in neutron matter and is given by the 
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expression 

deKK _8Gr,k /^VV • 2n 1^,2 f 2(l-7i)V[9(2 + n)] 
dco ~ Stt ^ l^mArJ ' 1 19/18 + II7I/9 + 27I/9 ' ^^^''^ 

where the upper hne corresponds to extra-dimensional gravitons, the lower line to dilatons, = / uj"^ 
and rrij = j"^/ is the effective mass of KK particles on the four dimensional brane; here Gat is Newton's 
constant, = +p^ and cos 9cm = p-p', where p and p' are the momenta of colliding neutrons, k and 
uj are the momentum and energy of radiated particles. The gravitational emissivity is obtained upon 
taking the phase-space integrals over the rates ()220|) and summation over the momenta j. The bounds 
obtained with the one-pion-exchange approximation and the full T-matrix (and an analog of the soft 
neutrino approximation) are 

f 7xl0-^mm (n = 2), 9 x lO'^mm (n = 3), , . 

\ 3 X 10-^mm (n = 2), 4 x lO'^mm (n = 3), ^ ^ 

where the upper and lower lines differ mainly in the treatment of the nuclear interaction: the upper line 
corresponds to the free space interaction between nucleon in terms of a T-matrix |221j . the lower line to 
the one-pion-exchange |222j interaction. These bounds provide one of the most important constraints 
on the size of extra dimensions, assuming that our current understanding of the supernova mechanism 
and energetics is correct. 



3.5 The role of pairing correlations in neutrino radiation rates 

Pairing correlations play a twofold role in the thermal evolution of neutron stars. At asymptotically 
low temperatures they suppress the neutrino emission processes exponentially, because the number of 
excitations vanishes as exp(— A/T), where A is the gap in the quasiparticle spectrum. At moderate tem- 
peratures T < Tc the pairing field lifts the constraints on the one-body (quasiparticle) bremsstrahlung 
(jl58|) and opens a new channel of neutrino radiation |1971 I198j . The corresponding charge- neutral 
current diagrams are shown in Fig. |27| These diagrams are associated with the following polarization 
tensor 

zn</^(g) = -2gB{qo) lmUv/A{q) = J2 J -^^A^" (P + Q)G^iP) TF<{p + q)F>\p)] , (222) 

where G{p) and F{p) are the normal and anomalous propagators, defined in Subsec. 12.2.21 the a- 
summation is over the spins. The superfluid in a neutron star can be considered as a two-component 
system, which, for a fixed density and temperature, consists of paired quasiparticles in the condensate 
and elementary excitations above the condensate. Their quasi-equilibrium densities are controlled by 
Cooper pair formation and pair breaking processes. The rate of these reactions are non-exponential 
in the vicinity of Tc, however, they are suppressed at asymptotically low temperatures exponentially, 
because of an exponential decrease of excitations above the condensate. These processes can proceed 
with the emission of neutrino pairs via the reactions {NN} — > N+N+u+u and N+N — > {NNy + u+u, 
where {A^A^} denotes the Cooper pair, A^ an excitation. Neutrinos of all three flavors can be emitted 
in such a process. The corresponding emissivities are 



G'pCNPFNm^TU ^ \ l[-l\e{TcN-T), N e {n,p) (223) 
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Figure 27: The one-loop contribution to the polarization tensor in the superfluid matter; 
solid lines refer to the baryon propagators, wavy lines to the (amputated) propagator. 



where C„ = ^17^ + ^25'i7^,o- and Cp = ^o7p + ^3- Here = 1, = if neutrons are paired in the ^Sq 
state and ^1 = 2/3 ^2 = 4/3 if they are paired in the state. Because of their relative low density 
protons pair in the ^5*0 channel and there is only a vector current contribution from protons; the factor 
^3 ~ 1 takes into account the electromagnetic correlations in the weak vertex. The functions 'jn/p and 
Injv^a take into account the correlations due to the strong force in the vector and axial vector vertices 
respectively. The temperature dependence of the emissivity eyy oc is familiar from the analysis 
of the reaction n— >-?T, + z/ + i/in magnetic fields [Eq. (jl96|) ] and is characteristic for the one-body 
bremsstrahlung process. The remaining parameters in Eq. (j223|) enter through the density of state of 
a single fermion. The integral in Eq. ()223|) is defined as 

/(a;) = r (""^^ ^)'^^ ^ a/^s-^^ (224) 

Jo [exp(x cosh y) + 1]^ V 4x 

It can be seen that in the limit A/T ^ 1, the rates of the CPBF process are exponentially suppressed, as 
it is the case for the competing two-nucleon processes. However, because of mild phase space restrictions 
(phase space volume of a single nucleon) these processes considerably contribute to the total neutrino 
emissivity at moderate temperatures T < T^. The magnitude of the CPBF processes is e^p ~ 

We now turn to the charged current weak decay Urea process n— i>p + e + z/in the superfluid matter 
and concentrate on the one-loop approximation |2UUj . This process is described by the first diagram 
in Fig. (jTTj) with the Z° replaced by and NN = n,p. The second diagram does not contribute at 
one-loop. 1 The vector and axial-vector one-loop polarization tensors have the form 

^n^M(?) = E / (0i^^(^ + ^)^^(^)' ^225) 
i.e. are identical for the vector and axial vector vertices; explicit evaluation of this expression leads to 

where = (1/2) (1 -|- ^p/ep), Up + = 1 with Sp = ^/^^ + and ^p being the proton single particle 
spectra in the superfluid and unpaired states; the quantities with the index k refer to the same functions 
for neutrons. Inspection of the denominators of four terms contributing in Eq. ()22(j|l shows that the 

^Note that ref. |2(JUj treats also the second diagram. While this diagram contributes at the second order in the strong 
interaction, and in the general case where the loops are summed up to all orders, it is strictly zero at one-loop order. 
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Figure 28: The neutrino emissivity in units of eo versus temperature (solid lines A„(0) = 
Ap(0) = 0.5 MeV, dashed-dotted lines A„(0) = 0.5, Ap(0) = 2 MeV). The scattering, pair- 
breaking contributions and their sum are shown by dashed and dashed-dotted and solid 
lines. The upper short-dashed line is the extrapolation of the rate for unpaired matter to 
low temperatures, the lower one corresponds to the exponential suppression as discussed in 
the text. 



first two terms correspond to excitations of a particle-hole pair while the last two to excitation of 
particle-particle and hole-hole pairs. The last term does not contribute to the neutrino radiation rate 
[uj > 0). We identify the first two terms as scattering {SC) terms, while the third term as pair-braking 
{PB) term. For unpaired neutrons and protons Up^k = 1 and Vp,k = 0, only the first term survives and 
one recovers the polarization tensor of non-superconducting matter. Upon evaluating the phase space 
integrals, the neutrino emissivity is written as 



eoJ, J 



) poo 

dy gsiy) [I''' + I''''] / dzz'Uz 

D Jo 



where eo is defined in Eq. p93p . This result differs from its normal state counterpart by the sum of 
the integrals /'^'^ + I^^ which reduces to the logarithm in Eq. ()193p for T > Tc. These integrals are 
defined as 



JSC 
jPB 



dx ( i""^ ^ ^^^/^^ 

up; 



x"^ + wl + y) 



(227) 



where x = /3^p, Wi = /3Ai {i e n,p), y = l3u and H-S'c/^^ = {u ± Ep) / ^/\(uJ'±~ep}^'^^'^^\; the exphcit 
form of the functions x^ is given in ref. |200j . Fig. |2H1 shows the temperature dependence of the direct 
Urea neutrino emissivity in the range 0.1 < T/T^ < 1. An important feature seen in this figure is the 
nearly linear dependence of the emissivity on the temperature in the range 0.1 < T/Tc < 1; the com- 
monly assumed exponential decay - a factor exp(— A/T) with A = max(A„, Ap) - underestimates the 
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emissivity. (Similar conclusion concerning the suppression of the direct Urea process by pair correlations 
were reached in ref. which treated the scattering contribution to the emissivity). The contribution 
of the pair-breaking processes becomes substantial in the low-temperature range 0.1 < T /T^ < 0.4 
where it is about half of the scattering contribution at T/Tc ~ 0.1. For unequal values of pairing gaps 
(e.g. Ap = 2, An = 0.5 MeV) the emissivity of the scattering process is suppressed, since at a given 
temperature the phase space accessible to the excited states is reduced. The pair-breaking processes 
are almost unaffected since they are related to the scattering of particles in and out-of the condensate. 

So far our discussion was confined to the case where the nucleons are paired in the ^5*0 channel. At 
high densities, as is well know, neutrons are paired in the channel. In this case the gap function is not 
isotropic in general and may have nodes on the Fermi- surface. The suppression of the phase space of the 
Urea process was studied in terms of suppression factors 71 = e^/eQ. In the case where the gap function 
has nodes at the Fermi-surface the suppression is a power law |151 1226j . Another possible implication 
of the P-wave superfluidity is the existence of Goldstone modes - collective excitations - which are due 
to breaking of the rotational invariance by the anisotropic P-wave condensate. These Goldstone boson 
or "angulons" (a) may couple to the weak neutral current |227j . The process a + a —>■ u + u leads 
to an emissivity which scales as eang oc 10^^Tg^(0.15/f )^ erg cm~^ s~^ where v is the angulon speed. 
The high power of temperature in the process is due to the thermal nature of angulons. At moderate 
temperature T < Tc the CPBF processes are more efficient than the angulon annihilation; at extreme 
low temperatures and in the cases where CPBF processes are suppressed exponentially, the angulon 
annihilation processes may become important. 

4 Cooling of neutron stars 

The equation of state of dense hadronic matter and the neutrino emissivities discussed in the previous 
chapters are the key ingredients of cooling simulations of neutron stars. The results of such simulations 
are combined with the experimental measurements of surface temperatures of neutron stars to gain 
information on properties of compact stars. In this chapter we review these simulations and their 
comparison with observational surface temperatures. 

Neutrino emissivities control the cooling rate of a neutron star during the first 10^ — 10^ yr of 
their evolution. For later times the photon emission from the surface dominates and the heating in 
the interior can be a significant factor in maintaining the surface temperature above the observational 
limit. Depending on the dominant neutrino emission process in the neutrino emission era t < 10^ yr 
the cooling proceeds according to the slow (standard) or the fast (nonstandard) scenario jTHl 1225] . 

The slow cooling scenario is based on neutrino cooling via the modified Urea and bremsstrahlung 
processes, modified appropriately to take into account the superfiuidity of the star's interior. The fast 
cooling scenarios invoke "exotic" cooling mechanism(s), such as the pion/kaon decay processes, the 
direct Urea process on nucleons or hyperons, as well as their counterparts in deconfined quark phase(s). 
We have seen that the phase-space arguments are the key to understanding the relative importance of 
different processes; those leading to fast cooling are characterized by a one body phase space and hence 
temperature dependences are ~ while those responsible for slow cooling originate from two-body 
processes and their emissivities scale with the temperature as T^. 

An inspection of the observational data on neutron star surface temperatures, which is commonly 
presented on a plot of photon-luminosity (or surface temperature) vs age (see Fig. 1221 and EOl below) 
shows that the data cannot be described by a single cooling track; there must be an effective mechanism 
that regulates the cooling rate in a manner that some of the stars cool faster than others. In particular 
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the Vela pulsar 0833-45 and the radio-silent compact star Geminga have temperatures that are far too 
low to be accommodated within a model based on slow cooling agents. It is reasonable to assume that 
the heavier stars cool via some fast mechanism, while the lighter stars cool slowly via the modified 
processes. The fast cooling agents are effective above a certain density threshold, which could be either 
the density of phase transitions to a novel state of matter or the density at which kinematical constraints 
are lifted, as is the case for the Urea process. Such an approach allows one to spread the cooling tracks 
within a range that can accommodate the currently available data. While such a strategy is a plausible 
ansatz one should keep in mind that neutron star cooling is a complex, multi-parameter problem and 
the processes which are actually responsible for fast (and to some extent slow) cooling of neutron stars 
are not known with certainty. 

The thermal evolution is governed by the coupled system of equations for energy balance \2'2H\ I229| 



- ,„v(-..e-^ke^-cf-^], (228) 



dr \ dt 

and thermal energy transport, 

9iTe^^') = (229) 
dr Airr'^K ' 

where L is the luminosity, T the temperature. This system requires as a microphysical input the 
neutrino emissivity e^ip, T), the heating rate h{p, T, Q, Q), with Q and Q being the spin frequency and 
its derivative, the heat capacity Cv{p,T), and the thermal conductivity, k{p,T), where p is the local 
density. The boundary conditions for Eqs. ()228j) and ()229|1 are L{r = 0) = and T(r = r^) = Te, where 
Te and Te are radius and temperature of the envelope. The last condition matches the temperature of 
the star to that of the envelope, since to a good approximation the temperature gradients and heat 
transport within the envelope are independent of the cooling of the "thermal core" of the star which 
extends from densities 10^° g cm~'^ to its center. After the first 10^ — 10^ years, the thermal core of the 
star is isothermal and the cooling of the star is described by Eq. (j228p which can be written in terms 
of quantities which are integrals of the microscopic parameters over the volume of the isothermal core 



jrri 

C,— = -L, -L., + H, (230) 
at 



where Lj, and are the neutrino and photon luminosities, H is the heating rate, and C^, is the 
integrated specific heat of the core. The photon luminosity is given by the black-body radiation formula 
Ly = AnR'^aT^, where a is the Stephan-Boltzmann constant, R is the radius of the star, is the surface 
temperature. The neutrino cooling era is characterized by the condition |Lj,| S> \L^\, \H\. The specific 
heat of the star is mainly due to degenerate fermions and scales with the temperature as C„ ~ T. If 
there is a single dominant neutrino emission process with known temperature dependence Eq. ()230p 
can be integrated; for typical slow processes (modified Urea, etc) L^, oc and one obtains t oc T~^; for 
rapid cooling processes (e.g. direct Urea) L^, oc and t oc T~^. The photon cooling era is characterized 
by the condition |Lj,| -C If the heating processes are ignored, the photon luminosity balances 

the thermal inertia term on the r. h. side of Eq. (|230|) . In more realistic cases the thermal inertia terms 
can be ignored and the photon luminosity is fully determined by the heating rate, which is a function 
of time. 



70 



4.1 Observational data 



Only a small fraction of radio pulsars that are visible through their radio emission have measurable 
photon fluxes from their surfaces or magnetospheres. A representative sample of 27 pulsars that will 
be used below to illustrate the constraints on the cooling theories by the observational data is given in 
Table ??. The effective surface temperatures are specified together with their 2a error range. 

The ages of the pulsars listed in Table ?? are their spin-down ages, r = P/2P, except the PSRs 
0531+21 (Crab), 0833-45 (Vela) and 0002-1-62, for which the age is know either through historical 
records (Crab pulsar) or the age of the supernova remnant they are embedded in. The spin-down age 
assumes that the star decelerates under the action of magnetic dipole radiation with a constant, time 
independent rate. If we write the star's spin-down rate as 

n{t) = K{t)n'^{t), (231) 

where Q is the spin frequency and n is the breaking index, the assumption above translates to K = const 
and n = 3. The spin-down age approximates the true age of a pulsar within an accuracy of a factor 
of 3 or so. The sample of pulsars in Table ?? is divided in four categories depending on a number of 
observational and fitting features. Four pulsars that belong to category D have not been detected in the 
soft X-ray range. The sensitivity of instruments sets an upper limit for the surface temperature. The 
data from ten pulsars of category C contain too few photons for spectral fits. The surface temperatures 
for these objects were obtained by using the totally detected photon flux. The spectra of eight pulsars 
of category B, which include the Crab pulsar 0531+21, can be fltted either (i) by a power-law spectrum 
or (ii) by a blackbody spectrum with high temperature and small effective area, much smaller than 
a neutron star surface. Presumably, their X-ray emission is dominated by magnetospheric emission. 
Therefore, the temperatures, determined from the spectral flts, are probably higher than the actual 
surface temperatures. Finally, there are flve pulsars of category A, 0833-45 (Vela), 0656+14, 0002+62, 
0630+18 (Geminga), and 1055-52, whose spectrum can be fltted by a two-component flt. The soft 
blackbody component is attributed to the surface emission, while the hard blackbody (or power-law) 
component is attributed to the magnetospheric emission. These pulsars are marked with error bars in 
Figs. EHl and Eni 

4.2 Cooling simulations 

Figure QUI shows the cooling of a family of neutron stars with masses in the range 1.0 to 1.9 Mq featuring 
the same microphysical ingredients P'3Dj . Baryonic matter is paired below the critical temperatures of 
the superfluid phase transition which are of the order of 10^ K; a phase transition to the pion condensed 
phase at the density ric = Suq is assumed. The Cooper pair-breaking/formation processes (CPBF) are 
not included in the simulation. The initial phase of the cooling for < t < 10^ yr is independent of 
the assumed value of the temperature at t = 0. At this stage, the star supports temperature gradients 
throughout the thermal core. The flrst kink in the cooling curves at t ~ 10^ marks the point where 
the cooling wave, which propagates outwards from the center of the star, reaches the surface. During 
the entire subsequent evolution the thermal core is isothermal and the slope of the cooling curves 
within the era 10^ < t < 10^ yr is determined by the dominant neutrino emission mechanism(s). 
The second kink at t ~ 10^ yr marks the point where the transition from the neutrino dominated 
cooling to photon dominated cooling occurs. The change in the slope is due to the difference in the 
temperature dependence of the photon and neutrino cooling rates. The neutrino dominated cooling 
era is independent of the initial temperature assumed at the instance when the star become isothermal 
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Table 5: Sample of observed data. The entries are: rotation period P and period derivative P, spin- 
down age T = P/'iP, K — PP and effective (redsfiifted) surface temperature T°^; bb and mH refer 
to blackbody and magnetized hydrogen atmosphere fits; f refers to known true age rather than the 
spin-down age. 
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Figure 29: Dependence of surface photon luminosity on the age [230]. The model includes 
pion condensation at high densities; the CPBF processes are switched off. The right panel 
includes pion mode softening effects in the modified Urea process |2()lj while the left panel 
is based on OPE result of ref. 



(unless this temperature is unrealistically low); in contrast, the initial conditions for the photon cooling 
dominated era strongly depends on the cooling pre-history of the star. The late time asymptotics of 
cooling in the photon dominated era is mainly determined by the models of the non-isothermal envelope 
which relate the temperature of the isothermal core to the surface temperature; the core temperature 
drops across the envelope by roughly two orders of magnitude. The late time asymptotics of cooling is 
likely to be dominated by the heating in the interior of the star. 

Let us now concentrate on the neutrino dominated cooling era. The qualitatively different behavior 
of stars with masses M > 1.6Mq and M < I.GMq arises from the fact that the former feature a pion 
condensate, since their central densities are larger than the threshold density for pion condensation. 
This segregation of the cooling tracks into high-temperature and low-temperature ones is common to 
the cooling theories with slow and fast cooling agents. An inspection of the observational data, shows 
that the neutrino dominated era cannot be fitted by a single track and a mechanism is needed to 
provide a smooth crossover from slow to rapid cooling. One possibility is the "threshold mechanism" 
in which some of the (low mass) stars have central densities below a threshold for a fast processes, and 
thus cool slow, while others, more massive, above this threshold and thus cool fast. Such mechanism is 
unsatisfactory, since one will need an extreme fine tuning of masses of neutron stars to accommodate the 
data (see e. g. Fig. 1221) • should be noted that any phase transition (as, e.g., pion or kaon condensation) 
or kinematically constrained process (such as the direct Urea process) will have precursors, where the 
fluctuations below the transition density will transform the sharp transition into a smooth crossover. 
The effects of such a precursor for the case of pion condensation, the so called pion mode softening, and 
its effect on cooling is studied in refs. |23()[ 1233] . 

A more interesting (and realistic) mechanism of the crossover from slow to fast cooling is offered 
by the CPBF processes. The variations in the cooling rates introduced by these processes arise due to 
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Figure 30: Dependence of surface photon luminosity on the age. The CPBF processes are 
included in all graphs. From left to right: the model includes pion condensation, the model 
includes in addition the Urea process; the model features neutron gap reduced by factor 

6. I23nil2SIl- 



the fact that the density profiles of pairing gaps for baryons map differently on the density profiles of 
light and massive stars. This effect is demonstrated in Fig. IHOl where the models presented in Fig. 1221 
are supplemented by the CPBF processes. The leftmost panel shows the effect of including the CPBF 
process in the presence of pion condensate, while the next plot to the right shows the effect of adding 
the direct Urea process. While the segregation between the high and low mass objects remain, the 
cooling tracks of low mass objects are now spread in a certain range. 

A wider range of surface temperatures can be covered upon suppressing (artificially) the neutron 
gaps by a factor of a few. In this case, the matter is unpaired at high-densities and the star cools as 

partially superfluid object. One may conclude that not only the magnitudes of the gaps are important, 

but also the density profiles over which they are spread [13 UHl l^32t I233j . 

Cooling simulations that do not include the heating processes in the interiors of compact stars fail to 
account for the highest temperatures in the sample of X-ray emitting pulsars. It is very likely that the 
late time evolution of compact stars including the entire photon cooling era and partially the neutrino 
cooling era are significantly affected by the heating processes in the star's interior. The heating processes 
can be roughly divided into three categories: (i) the heating due to the frictional motion of neutron 
and proton vortices in the superfluids; (ii) the heating that arises due to the local deviations of matter 
composition from /3-equilibrium; (iii) the heating due to mechanical processes in the solid regions of the 
star (crust cracking, etc.). Describing these processes will require a detailed account of the physics of 
neutron star interiors on the mesoscopic scales and is beyond the scope of this article (see ref. |231j for 
details). 
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5 Concluding remarks 

This review covered a number of topics that are relevant to the theory of compact stars, in particular the 
physics of hadronic matter with baryonic degrees of freedom and weak interaction processes involving 
hadrons. It should be clear from the presentation that despite the enormous progress achieved during 
the past decade, the theory is far from being completed and a large number of exciting topics remain 
to be studied in the future. 

Compact stars continue to pose an enormous intellectual challenge to the physics and astrophysics 
communities. Rapid progress during the last four decades since the discovery of the first pulsar is 
an evidence of the vitality of this field. Current and planned observational programs continue the 
exploration of compact objects in the electromagnetic spectrum; studies of compact stars through 
gravity waves may become possible in the near future. All this provides an excellent basis for future 
theoretical studies of compact stars. 
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